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CHAPTER I 


A VECTOR AUTOPILOT SYSTEM 


1-1 INTRODUCTION 

An essential requirement of an aircraft attitude con- 
trol system is that deviation of the body axes relative to 
a reference axes frame must be sensed. In addition, to 
overcome the ever-present possibility of errors or failure 
of the sensors, various configurations of redundant sensors 
are usually employed to assist in detection and correction 
of errors. To this end, there has been a continuing effort 
to improve existing sensors, to develop new sensor config- 
urations, and to develop new sensor devices. 

This chapter discusses the role of a vector magneto- 
meter^ as a new instrument for aircraft attitude determin- 
ation. Although magnetometers have played a role in the 
attitude measurement of missiles' and satellites [Ref. 1-1], 
there is an apparent lack of application in aircraft systems 
By providing independent measures of attitude, the solid 
state vector magetometer sensor system can not only assist 
in improving accuracy and reliability of existing systems 
but can also reduce component count with obvious benefits 
in weight and cost. Additionally, since a large number of 
aircraft heading reference systems depend on measurement 
of the Earth’s magnetic field, it can be shown that by sub- 
stituting a three-axis magnetometer for the remote sensing 
unit; both heading and attitude measurement functions can 
be derived using common elements, thereby further reducing 
the component count. 


^Aviation use to date has been essentially scalar magneto- 
metry . 
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To investigate the feasibility of the above system, this 
chapter will proceed by developing a technique to determine 
attitude given magnetic field components. Sample calcula- 
tions are then made using the Earth's magnetic field data 
acquired during actual flight conditions. Results of these 
calculations are compared graphically with measured attitude 
data acquired simultaneously with the magnetic data. The 
role and possible implementation of various reference angles 
are discussed along with other pertinent considerations. 
Finally, it is concluded that the Earth's magnetic field as 
measured by modern vector magnetometers can play a signifi- 
cant role in attitude control systems. 

1-2 ATTITUDE DETERMINATION 


Coordinate systems are usually defined by orthogonal 
right-handed sets of three unit vectors. An example of such 
a set is illustrated in Fig. 1-1 where the orientation of 
the body fixed frame used in this paper is delineated. An- 
gular' rotations are conventionally defined as rotations in 
the plane normal to a unit vector with the positive sense of 
rotation defined by the right-hand rule [Ref. 1-2], 

To derive relationships of attitude variations as a 
function of magnetic vector component variation, we can pro- 
ceed by considering matrix representations of an orthogonal 
transformation. If Hx, Hy, and Hz are the magnetic compon- 
ents measured at a desired airframe attitude and Hx', Hy ' , 
and Hz' are the components measured after any rotation of 

the body, vector H' = [Hx’ Hy' Hz ' ] can be related to 

T 

vector H = [Hx Hy Hz] by an orthogonal linear transform- 
ation H' = AH. Here A must satisfy the orthogonality con- 
T T 

dition AA =1, where A is the transpose of A; addition- 
ally, the determinant of A must be unity [Ref. 1-3, 1-4] . 
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Fig. 1-1 AXIS ORIENTATION 

Rotations about the z axis in Fig, 1-1 result in yaw 
deviations (^) and in new components (H'), as shown by 


Hx' 


cos sin fj; 0 
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Similarly, independent rotations about the y axis and 
the X axis result in pitch (0) and roll C'^) dependent varia- 
tions in the measured H components, as shown by 
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The effect of a combined rotation can be expressed by 
using the product of the transformation matrices. In add- 
ition, if the rotations are small, the total rotation exper- 
ienced by applying sequential rotations is independent of 
the order in which the rotations are performed [Ref. 1-3,1-4J. 
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Assume that the angular variations 0 , ip , and <j) are 
small enough so that the small angle approximations 

sin 0 S' 0 , sin ^ , sin , 

cos 0 == cos i|i cos <j) S' 1 

can be made. Then, if the products of small angles (in rad- 

ians) can be assumed to be much smaller than the angles 
alone, the expression reduces to 

"Hx' 1 r 1 i(< -el THx" 

Hy' = -ti) 1 <p Hy 

_H2'J L 9 ij Lh2:J (1-5) 

Further modifications in the form of the matrices result in 

"Hx'" '-Hz HyOl fSl PHx' 

Hy' = 0 -Hx Hz i|j + Hy 

.Hz'J L Hx 0 -HyJ [_ <^ J L^ 2 j (1-6) 

By subtracting, . we arrive at an expression for the dif- 

ference in H components as functions of angular deviation. 

'Hx' 1 ~AHx' f-Hz Hy 01 fS* 

Hy' - Hy = AHy = 0 -Hx Hz ijj 

_Hz'J LhzJ L^HzJ [ Hx 0 -HyJ L-j)J (1-7) 

It is significant to note at this point that the trans- 
formation matrix is singular implying that solutions for 
9, and'(i) are not independently available. 
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1-3 ATTITUDE DETERMINATION EMPLOYING MAGNETIC 
FIELD COMPONENTS 


A given orthogonal set of three unit vectors can be dis- 
placed in Euclidean space by rotating the system through any 
angle 6 about a directed rotation axis. It is also customary 
to represent this rotation vectorially as a directed line 
segment whose length is proportional to the rotation angle. 
This rotation is analogous to the rotation experienced by the 
body fixed frame of Fig. 1-1 as the aircraft experiences 
combined pitch, yaw, and roll variation. During flight the 
body fixed set rotates about this rotation axis assuming new 
(possibly erroneous) attitudes in space. The task of the 
attitude sensing system is to provide measures of compounded 
pitch, yaw, and roll that would result in the same attitude 
assuming that the rotations occurred sequentially about the 
X, y and z axes rather than the actual rotation axis. 

It was shown in the previous section that a compounded 
rotation of an orthogonal set can be described by a product 
of respective transformation matrices. Additionally it was 
noted that for small angular rotations the order of multi- 
plication is unimportant. Using the relationships of (1-7) , 
expressions for the angular deviations in terms of measured 
magnetic vector components can be derived. 

AHx == -Hz8 + (l-8a) 


yields . 


6 = (Hyi|» - AHx) /Hz 


(l-8b) 


^ = (AHx + Hz0)/Hy 


(l-8c) 
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Similarly, 

AHy = -Hxi|; + Hz()) Cl-9a) 

yields 

ij/ = (Hz<^ - AHy}/Hx 

a = (AHy + Hxi(j)/Hz 

and 

AHz = Hxe - ny<^ (1-lOa) 

yields 

e = (AHz + Hy 4 .)/Hx (1-lOb) 

= (Hx3 - AHz)/Hy (1-lOc) 


(l-9b) 

(l-9c) 


Assuming that Hx, Hy and Hz are nominal vector compon- 
ents as measured in a reference attitude and that Hx ' , Hy ' 
and Hz' are new field components at the new attitude, then 
AHx = Hx ' - Hx, AHy = Hy' - Hy, AHz = Hz' - Hz are expres- 
sions of the incremental changes in field components. Add- 
itionally, before using (1-8) , (1-9) or (1-10) to solve for 

attitude variations (pitch, yaw, or roll) , one additional 
angle from an auxiliary sensor^ must be supplied. Using one 
additional angle of rotation (about any one axis) the remain- 
ing two rotations can then be calculated. 

To illustrate this point, flight data acquired during 
the flight of a NASA flown Convair 900 instrumented with a 
three-axis magnetometer and a Litton inertial navigation 
system were used to calculate roll, pitch, and yaw. 


^It was noted following (1-7) that a unique solution for 
attitude variation is not possible using magnetic field 
data alone. 



Attitude variation about each of the three axes was calcul- 
ated using measured magnetic field components supported by 
one angle from the inertial system. The results of these 
calculations are plotted in Figs. 1-2 through 1-4. 

It is significant to note that the rotations shown 
occurred simultaneously (i.e./ time base is the same for all 
three figures) . The flight was at an altitude of approxi- 
mately 5000 ft at an airspeed of approximately 250 nmi/h. 

Although the data used to plot the attitudes shown in 
Figs. 1-2 through 1-4 were not acquired specifically for this 
purpose, the correlations in measured and calculated atti- 
tude clearly show that*, within the limits of instrument ac- 
curacy, signals proportional to attitude variation can be 
derived using flight data. 

1-4 A POSSIBLE SYSTEM CONFIGURATION 


•Since the intent of this chapter is to introduce the 
notion that magnetometer technology has advanced to the point 
where three-axis magnetometers can be incorporated in air- 
craft attitude sensing systems on a cost effective basis, the 
system discussion will be limited in scope to describing a 
possible combined heading and attitude measurement method. 

Heading references fall into three classes; 1) those 
that depend on the Earth’s magnetic field, 2) those that de- 
pend on the use of low-drift gyroscope to retain a preset 
azimuth, and 3) those (gyrocompasses) that depend on sensing 
the Earth's rotation [Ref. 1-5]. By far the greatest number of 
aircraft heading systems depend on the Earth's magnetic 
field, although many of these include gyroscopes to' improve 
the performance characteristics. 
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FIGURE 1-3 PITCH AXIS 




FIGURE 1-4 YAW AXIS 
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A popular system combination (with no gyro) is to 
combine a pendulous remote magnetic sensor and a synchro 
receiver in a null seeking circuit. The philosophy being to 
attempt to measure only the horizontal component of the 
Earth's magnetic field and to swing the receiver into align- 
ment with it. Under acceleration, departures of the sensor 
unit from the horizontal result in angular heading errors e 
[Ref. 1-5]. 

£ = (aH/g) tany sine 

where aH is the horizontal acceleration, g is the accelera- 
tion due to gravity, 9 is the angle between the acceleration 
vector and magnetic north, and y is the magnetic field dip 
angle; arctan (vertical field/horizontal field). 

Accuracy of this system can be improved by incorporat- 
ing a strapped-down solid state magnetic sensing unit (free 
of acceleration errors) that measures and displays the angle 
of the Earth’s horizontal magnetic component relative to the 
aircraft. This system can be implemented as follows: 

1) Determine the direction of the magnetic vector F rela- 
tive to the sensors (and the airframe), by measuring the x, 
y and 2 components (Figs. I-l and 1-5). The direction co- 
sines cosa, coss, cosy are the cosines of the angles a, S, y 
between the magnetic vector and the positive x, y and z 
axes. Additionally, 

.1 

cosa = x/(x^ + t 2 ^)'^ 

1 

cosg - y/(x^ t t z^)^ 


1 


cosy = z/(x2 + y2 + z^)^ 
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2) Using either a vertical reference^ or knowledge of 
aircraft attitude / we can effectively rotate the body axes 
such that the x-y plane is horizontal (see Chapter II) . 

3) Simple application of direction cosines will yield the 
direction of magnetic north in the aircraft's x-y plane. 


X 


z 


Fig. 1-5 FIELD VECTORS AND DIRECTION COSINES 

Although the preceding discussion implies that heading 
can be determined by using a strapped-down magnetometer, 
there remains the problem of attitude determination. An- 
other widely used system for obtaining a heading reference 


^Not necessarily -derived inertially [Ref. 1-11]. 
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is to combine the relatively excellent short term stability 
of a directional gyroscope with the long term stability of 
magnetic field measurements. By slaving the directional gyro- 
scope to the magnetic field [Ref. 1-5, sec. 10.4.7], gyroscopes 
with relatively large free drift error can be used to pro- 
vide an excellent heading reference. 

Replacement of the pendulous remote sensing unit of this 
type of system with a strapped-down vector magnetometer would 
result in both heading and attitude information on a contin- 
uous basis. This combination would operate as follows; 

1) The system is initialized by determining a reference 
attitude (perhaps by using a primary inertial attitude sys- 
tem) . 

2) The angular position of the horizontal magnetic field 
component is computed as above and used to slave the direct- 
ional gyroscope. 

3) The directional gyroscope, with relatively good short 
term stability (devices with free drift of less than 0.5 
deg/h have been designed) , is used to determine yaw (>p) 
errors . 

4) For small angle deviations, (1-8), and (1-9), and 
(1-10) can be employed to recalculate aircraft attitude. 

The process loops back to step 2) closing the loop on a 
combined attitude and heading reference system. 

The sampling frequency required to maintain an accept- 
able level of error is of course determined by the aircraft 
performance expected (angular rates) and by the gyro error 
(drift rate plus errors due to additional sources such as 
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gyroscope tilt from vertical) . The overall system is such 
that heading can be determined as before with errors due to 
sensor departures from hori 2 ontal substituted for long term 
accumulation of attitude uncertainty (this can be corrected . 
by looping to step 1) at a frequency dependent on error 
rates) . Additionally one gains measurements of attitude 
with minimal computation and replacement of a mechanical re- 
mote sensing unit with a solid state strapped-down magneto- 
meter sensor. 

1- 5 OTHER CONSIDERATIONS 


The characteristics of the Earth's magnetic field and its 
variations have long been established [Ref. 1-6-1-10]. Since 
the field is to be used as a reference in the attitude 
measurement scheme, there is a need here to discuss its ad- 
verse characteristics.- Although the field does experience 
variation, most of the variation is either in amplitude 
(ionospheric contributions) or has time constants that make 
the variation negligible (secular variation) . 

In traversing local anomalies, there will, however, be 
deflections in the ambient field due to the additive effect 
of local dipoles or monopoles. The effect of local terrain 
caused anomalies can be visualized by picturing the main 
field vector oriented in space with a second modulating 
vector rotating at its tip. Maximum angular error would 
occur when this modulating vector has maximum magnitude and 
is positioned at right angles to the main vector. 

To illustrate the effect of local anomalies one can cal- 
culate the level of anomaly required to cause an error. 

Since the Earth's main field is typically in the order of 
0.50 G it is readily apparent that a local anomaly of ap- 
proximately 0.01 G at right angles to the local field is 
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required to cause an error of 1 deg. Furthermore, the 
local anomaly would have to be aligned with one of the air- 
craft body axes to result in one degree of attitude error in 
any one axis. Fortunately, anomalies with components of 
this magnitude positioned at right angles to the main field 
are extremely rare. In addition, the anomalies are local- 
ized over ore bodies or other geophysical irregularities, 
have magnitudes that diminish as the cube of altitude, and 
tend to average to zero over relatively short distances. In 
summary, the probability of encountering an anomaly that 
would cause as much as a 1 degree error is relatively small. 
The error, if introduced, will be short lived and, unlike 
drift error, will average to zero. 

Fundamental to a magnetic field referenced system is 
the ability to measure orthogonal components of the field 
vector. Precision and accuracy of measurement of the com- 
ponents is of course specified by the desired control spec- 
ifications . 

Since the Earth's magnetic field varies in magnitude on 
a global basis between 0.3 G and 0.6 G (30,000 gamma to 
60,000 gamma), it is apparent that full scale measurements 
of 0,6 G can be expected. Sensors mounted at right angles 
to the field will monitor no measureable field and thus de- 
fine the lower limit of measurement to be zero. For the con- 
tinental United States the declination varies between 60 and 
80 deg, resulting in a range in horizontal component of 0.15 
to 0.25 G with vertical component in the range of 0.4 to 0.55 
G. Heading variations (yaw) result in changes of the hori- 
zontally sensed field components and would specify the maxi- 
mum precision required. In addition, flight at 45 deg ± (n 
X 90 deg) (where n is any whole number) with respect to mag- 
netic north results in minimum sensitivity of the x and y 


17 


axes measurements. In this case sensor inputs would range 
between 0.106 and 0.177 G with minimum field at the north. 
Assuming the preceding ambient measurements, variations in 
component magnitude of approximately 0.0180 to 0.0305 G/deg 
for small angle variations can be expected. 

A brief survey of commercial magnetometer manufacturers 
reveals that triaxiai magnetometers that measure from zero 
to 0.6 G with linearities of 0.5 percent, noise less than ±1 
mG and sensitivies of at least 2.5 V per 600 mG are current- 
ly available. In addition, these devices have a bandwidth 
of direct current to at least 500 Hz and are rated to have 
less than 1 deg error in orthogonality. 

From a precision standpoint, it is apparent that vari- 
ations in yaw for this worst case situation can be sensed to 
better than 0.1 deg with currently available magnetometer 
technology. The sensor technology required to implement an 
attitude sensing system of reasonable specifications is avail- 
able '(more detailed analysis is presented in Chapter III). 

Although the preceding calculations indicate that for 
small angular variations attitude can be calculated using 
measured magnetic data, there is a need to consider the 
effects of larger finite rotations. In this case the small 
angle assumptions would not be valid and an Euler transform- 
ation would have to be made. Measurement of three axes of 
field components could be used to develop the direction co- 
sines required to determine the orientation of the axis of 
rotation, the angular rotation about it, and the three .angu- 
lar rotations of pitch, roll, and yaw. 

For the special case where the axis of rotation aligns 
with the magnetic vector, there would of course be no 
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measured component changes.'* By measuring the attitude of a 
second vector (not in alignment with the magnetic vector) , 
we could resolve the ambiguous situation cited above and pro- 
vide additional redundancy. 

The optimum auxiliary vector would be one that could be 
sensed without using inertial devices. The Earth's electric 
field can be considered. The main reason for considering 
this field as a means of providing an auxiliary angular ref-' 
erence is that the resultant system has the potential of be- 
ing completely solid state. The electric field vector can 
be used to determine attitude variation in a manner analogous 
to the magnetic vector system. Inherent limitations of each 
single vector system can be obviated if the vectors are not 
coincident. 

Although Hill [Ref, 1-11] reported success in controlling 
pitch and roll using the electrostatic field alone, comments 
by Markson [Ref. 1-12]. indicate that the electrostatic field 
is not always a reliable vertical reference. Employment of 
the electrostatic field for this attitude measurement system 
is limited to augmenting the magnetic field measurements by 
eliminating ambiguity of motion around the magnetic vector. 

The requirement of vertical electrostatic field is thus remov- 
ed and replaced by a requirement that the field direction is 
relatively stable. 

By using two independently derived vectors we have suf- 
ficient data to obviate the ambiguity just cited and we have 
the potential of providing redundancy as well. 


^*An example of this would be yaw rotation while flying 
straight and level over the magnetic poles or roll rot- 
ation while flying towards a pole at the magnetic equator. 
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1-6 CONCLaSION 

This chapter has identified a novel method of measuring 
aircraft attitude using relatively inexpensive, well devel- 
oped instrumentation. It has recognized that magnetic field 
sensing systems have been used to some extent in attitude 
sensing and control of space vehicles; it has also suggested, 
however, that with appropriate support, magnetometers can 
find increased application in aircraft attitude measurement 
systems . 

This claim is corroborated by actual flight test data. 
Magnetometers have evolved to a point where three axis mea- 
surements of the Earth's magnetic field can be made with 
sufficient precision and accuracy to enable measurement of 
small angle attitude variations. 

This chapter has also discussed a possible system con- 
figuration combining heading determination and attitude 
measurement functions. By replacing the conventional re- 
mote sensing unit with a three- axis magnetometer, it has 
been suggested that both functions can be obtained' with the 
hardware required previously for heading measurement alone. 

As with any system, there are limitations imposed. The 
main limitation for a vector magnetometer system seems to be 
the inability to sense rotations around the magnetic vector 
itself. This problem is not unlike the ambiguity experienc- 
ed by magnetic heading systems at high latitudes. By judic- 
iously incorporating auxiliary instruments, not only can 
the ambiguities be removed but a degree of redundancy can be 
added while still maintaining a cost and weight advantage 
over comparable systems. 


CHAPTER II 


AN ATTITUDE INDEPENDENT REMOTE MAGNETIC INDICATOR 
2-1 INTRODUCTION 

Preliminary investigation [Ref. 2-1] revealed that air- 
craft attitude can be calculated using measurements of earth's 
magnetic field vector and a single auxiliary rotation angle. 

An algorithm to compute the two remaining aircraft rotational 
angles was developed. Using flight data, it was demonstrated 
that an excellent correlation in computed versus actual air- 
craft attitude could be achieved. In addition to providing 
measurements of the magnetic field for redundant attitude 
computations (to improve accuracy and reliability of existing 
autopilot systems) , it was noted that the vector magnetometer 
could substitute for the remote magnetic sensing unit. In 
this manner both heading and attitude measurements could be 
derived using common elements with obvious benefits in weight 
and cost. 

This chapter discusses the mechanization of a micropro- 
cessor based computer system that uses a three axis magneto- 
meter plus gyro data to compute heading. The magnetometer is 
a three axis solid state device that can be mounted in a 
strapped down configuration resulting in an attitude independ- 
ent remote magnetic indicator. Gyro measurements of pitch and 
roll angle plus three axis magnetic measurements are used by 
the algorithm to compute aircraft heading. The system can 
function independently to compute heading or by simply in- 
creasing the stored program could implement the attitude com- 
puting algorithm of [Ref. 2-1] as well. 

The chapter proceeds by developing an algorithm to com- 
pute aircraft heading using the strapped down magnetometer 
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and two gyro measured angles. Practical aspects of designing 
the system including both hardware and software are then pre- 
sented. In addition, the limitations in instrument accuracy 
and operation as determined by sensor errors, signal process- 
ing errors, arithmetic precision and computation speed are 
discussed. Considerable computational capability inherent in 
the system enables minimization of systematic errors. It is 
demonstrated that inexpensive sensors can be employed with 
offset and orthogonality errors compensated by microprocessor 
programming. Finally, it is concluded that a microprocessor 
based computer with a solid state magnetometer can play a sig- 
nificant role in aircraft instrumentation. 

2-2 AN ALGORITHM TO COMPUTE AIRCRAFT HEADING 


Coordinate frames are usually defined by orthogonal 
right-hand sets of three unit vectors. An example of such a 
set is illustrated in Fig. 1-1 where the orientation of the 
body fixed frame used in this chapter is delineated. The ref- 
erence coordinate frame referred to in this chapter is orient- 
ed with axes x and y in the horizontal plane and axis z verti- 
cal (z down is positive). Pitch attitude angle (0) of an air- 
craft is defined [Ref. 2-2] as the angle between some prefer- 
red longitudinal axis and the horizontal reference. In this 
chapter, pitch angle is the angle between the x axis of the 
aircraft and the x-y plane of the reference axis set. Since 
angular rotations are conventionally defined as rotations in 
the plane normal to a unit vector with the positive sense of 
rotation defined by the right-hand rule [Ref. 1-2] , we will 
define positive pitch angle (0) as the ’’nose up" or positive 
rotation about the y axis when the y axis is horizontal. The 
roll and yaw angles ($ and f) will then simply be rotations 
about the x and y axes respectively. 
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By aligning the three magnetometer axes with the respec- 
tive X, y and z axes of the aircraft, we can measure magnetic 
field components of the aircraft at any attitude. For the 
trivial case where pitch (0) and roll ($) are both zero de- 
grees, Hx and Hy are the horizontal field components and we 
can compute yaw from the horizontal vectors as follows: 

H'l = cos Mhx/(Hx^ + Hy^)^ (2-la) 


or 




sin“^ (Hy/(Hx2 


Hy2) ^ 


(2-lb) 


We select either (2-la) or (2-lb) based on the relative mag- 
nitudes of Hx and Hy. By minimizing the numerator of the 
argument we guarantee that the inverse trigonometric operation 
results in an angle between zero and forty-five degrees with 
maximum sensitivity ensured. Heading is then computed using 
the signs of Hx and Hy to select the appropriate equation from 
Table 2-1. 



Table 2-1. Formulae to Compute Heading 

For most cases, the pitch and roll angles are not zero 
and inverse rotations are required to determine the actual 
horizontal field components Hx and Hy. Since any aircraft 
attitude can be represented as a sequence of rotations about 
each axis beginning at some reference attitude, we can 
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• . the reference Hx and Hy field components by perform- 

roll computation can be developed by cons id- 
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i, = X2(i2 • it) + yt (j2 • it) - 

A 

H • i2 = X2 


H 


(2-4a) 

(2-4b) 


and from (2-3) we obtain 


H 


• i9 = X3 (is 


i^ * io) ^33 * ^2 


L) + 23(^3 • ^ 2 ) 


H • i2 = ^3 


(2-5a) 

(2-5b) 
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then 


Similarly, 


X 2 = X 3 


(2-6) 


and 


H* j 2 = Y 2 = X3U3 . 32) + 73(33 • 32) + 23(^3 • 32) 

(2-7a) 

Y2 “ Y3 cos <5-23 sin 

(2-7b) 

H • ka = 22 = X3 U3 • i<2) + Y3 ( j 3 • ^2) + Z3 (1^3 • ka) 

(2-8a) 

Z 2 = Y 3 sin <i' + 23 cos 


These expressions can be summarized as 


{2-8b) 
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Fig. 2-1 AXES ROTATED IN ROLL 
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Similarly, considering an axis set rotated in pitch as shown 
in Fig. 2-2, we can express the reference set x^ , yj , in 
terms of the rotated set X 2 / Y 2 , ^2 follows 




cos 0 0 sin 0 ' 


^X2" 

yi 


0 10 

• 

Yz 

.^1 , 


-sin 0 0 cos 0 _ 


_ ^2 _ 



Fig. 2-2 AXES ROTATED IN PITCH 

Finally, if we assume that the axis set subscripted with ^ 
represents components of Earth's magnetic vector measured at 
an arbitrary aircraft attitude, we can derive the magnetic 
components (Hxh, Hyh, Hzh) in the horizontal plane for a 
given heading 



~ HKh" 


cos 0 0 sin 0 ‘ 


"10 0 " 


“HX3" 

Hyh 

= 

0 10 


0 cos $ -sin $ 

• 

HY3 

_Hzh _ 


_-sin 0 0 cos 0_ 


_0 sin $ cos 


1 — 

N 

to 

1 


(2-lla) 


“Hxh' 


cos 0 (sin 0 sin $) sin 0 cos 


'HX3- 

Hyh 

= 

0 cos $ -sin $ 

• 

Hys 

_Hzh_ 


j-sin 0 (cos 0 sin $) cos # cos 9_ 


-Hz3_ 


(2-llb) 


The algorithm to be implemented with the microprocessor 
would therefore require operations as outlined in Fig. 2-3. 
Details of programming method/ modifications to the above 
equations to facilitate programming and computation speed 
versus accuracy tradeoffs are discussed in following sections. 


DO FOREVER 


Measure, Digitize and Store 
Hx, Hy, HZ/ 0 and 


Correct Sensor Errors 


Compute the Horizontal Field 
Components Using Equation 2-11 


Compute Heading Using Equation 2-1 


Display Heading 


Fig. 2-3 LOGICAL OPERATIONS REQUIRED TO COMPUTE HEADING 
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2-3 MECHANIZATION OF THE HEADING ALGORITHM 

A. General Considerations 

To evaluate the performance of an integrated system ex- 
perimentally, an instrument was designed to implement the al- 
gorithm developed above. Several approaches were considered 
to implement the heading instrument for experimentation: 

1} A minicomputer implementation incorporating an 
HP-2100 minicomputer supported by peripheral interface and an- 
alog circuitry. Programming of the HP-2100 would have enabled 
the computer to. control multiplexing and processing of sensor 
data as suggested by Parish and Lee [Ref. 2-3] . 

2) A hybrid system composed of a remote data acquisition 
system to collect data from sensors for subsequent processing 
by a computer (possibly an HP-2100) . 

3) A digital/analog electronic implementation incorpor- 
ating the design of a special purpose computer to perform the 
required functions of a heading instrument. 

The first two approaches were abandoned since it was de- 
sirable to perform the experiments at various locations remote 
from a computer facility and to have data available immediately 
without having to rely on off-line computations at a later date. 
The design task then evolved to the design of a special pur- 
pose computer system to implement the algorithm, provide a 
means for evaluating the performance of the proposed algorithm 
and to allow modifications to the system if required. 

B. Design Criteria 


Having decided on the general approach to implementing 
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the algorithm it became necessary to consider the performance 
criteria desired of the instriment. 

1) Accuracy 

As a design goal, an absolute accuracy of +1.0° in 
heading uncertainty was selected for the laboratory 
implementation. This accuracy is compatible with 
coramercially available heading systems. 

2) Computation Speed 

The bandwidth of the system is determined mainly by 
the computation speed of the computer^. As a design 
goal, complete heading updates once per second was 
established. 

3) Flexibility 

A desirable feature of the laboratory evaluation 
instrument was considered to be flexibility. Re- 
visions or additions to the algorithm as predicted 
by experimental data should be incorporated with 
minimal redesign of the instrument. 

2-4 CONCLUSIONS 


An instrument designed to implement the heading algorithm 
developed above uses a three axis magnetometer to measure mag- 
netic field data in the vicinity of an aircraft. Since the 
magnetometer proposed is a solid state three axis fluxgate de- 
vice and is permanently mounted in a strapped down configura- 
tion, the implementation results in an attitude independent 


^The response times of the various sensors and analog circuitry 
are orders of magnitude greater than the desired one second 
sample interval. 
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remote magnetic indicator^. 

Several factors will contribute to system inaccuracy. A1 
though the major error sources can be evaluated mathematically 
j(Chapter IV) , there is a need to evaluate the implementation 
experimentally. Systematic errors that arise can be reduced 
by instrument computation. This capability (inherent with a 
computer based system) enables incorporation of less expensive 
sensors in the heading instrument with less concern with fact- 
ors such as temperature regulation, sensor orthogonality and 
sensor offset‘s. 

Since the algorithm can be implemented using a micropro- 
cessor as the major computer element, the resulting instrument 
will have inherent computation capability, be small in size 
and consume relatively little power. These factors, make the 
instrument an ideal device for aircraft application where the 
need for redundant distributed processing capability is in- 
valuable . 


^Current remote magnetic indicators are pendulous and rely on 
gravity to enable measurements of the horizontal magnetic 
vector (not attitude independent) . 

^Assuming that the sensors have repeatable or measurable char- 
acteristics, algorithms can be developed to correct previous- 
ly measured erroneous data. 
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DESIGN OF A MICROPROCESSOR BASED HEADING INSTRUMENT 
3-1 INTRODUCTION 


Progress in device and component technologies during the 
1970 's has led to an assortment of sophisticated integrated 
circuits (IC) devices [Ref. 3-1] which enable the design of in- 
struments with a high degree of sophistication and accuracy. 

Of these devices, the microprocessor has to date been the most 
exploited component in industrial control and instrumentation 
applications [Ref, 3-2 through 3-7]. There have been many 
papers presented addressing the general application and feas- 
ibility 'of applying microcomputers to particular design tasks 
[Ref. 3-8 through 3-21]. 

Although much of the literature to date on microproces- 
sors has addressed the design of commercial products (usually 
the final result of a carefully orchestrated effort beginning 
with a market survey) , the design of a laboratory instrument 
for algorithm evaluation differs in design philosophy. In 
particular, the laboratory instrument is designed to evaluate 
a proposed algorithm under laboratory conditions. The tradi- 
tional benchmark evaluations and attempts to match the micro- 
processor to the application is not only difficult but unnec- 
essary, If the processor is much more powerful than necessary, 
the "overkill" is little noticed; but if an insufficiently en- 
dowed microprocessor is selected, the effects can be devastat- 
ing. Not only will the program be difficult to write and vor- 
acious of memory, it would be difficult to change to a more 
powerful microprocessor part way through the project. With 
these considerations in mind, a general purpose, flexible 
microprocessor with powerful architecture and instruction set 
the Signetics 2650 microprocessor -[Ref. 3-22] was selected. 
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3-2 HARDWARE DESIGN CONSIDERATIONS 


The design of a microprocessor based system begins by 
considering the total system level block diagram to be imple- 
mented (Fig. 3-1). Inputs from five sensors including x, y 
and 2 axis magnetic data plus pitch and roll angles (Hx, Hy, 
Hz, 0 and <l>) are to be multiplexed, sequentially sampled and 
converted to a digital representation prior to processing (ex- 
ecuting the algorithm developed above). The main subsystem of 
Fig. 3-1, the central processing unit (CPU), operates under 
control of instructions stored in the system memory and inter- 
faces with the input and output subsystems via data ports. 

At this early stage in the design, it is significant to 
note that the block diagram of Fig. 3-1 differs slightly from 
that of a classical discrete hardware solution. The input 
subsystem (composed of analog multiplexer, sample and hold, 
and analog to digital converter) differs from a conventional 
data acquisition in that it is devoid of a control section. 

The microprocessor will control the data acquisition sampling 
and conversion in addition to performing the arithmetic func- 
tion associated with the algorithm. 

Having established a tentative block diagram of the in- 
strument, the design continues by addressing relevant charac- 
teristics and limitations of each subsystem. These character- 
istics will then in term be considered in configuring the 
final system and program to be executed. 

1) The Analog Subsystem 

Composed of the analog multiplexer, sample/hold and an- 
alog to digital converter, the analog subsystem of Fig. 3-1 
affects both system accuracy and throughput rate. The 


THE ANALOG SUBSYSTEM 
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Fig. 3-1. SYSTEM BLOCK DIAGRAM 
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well-known Shannon theorem [Ref. 3-23, 3-24] on sampling 
theory defines one of the basic limits on throughput rate 
stating that the minimum frequency for sampling must be double 
the highest significant frequency of the signal, including the 
noise on the signal. This minimum frequency is necessary, the 
theorem states, if the sampled signal is to contain all of the 
information needed for undistorted reconstruction. At a lower 
sampling frequency ailiasing can occur^ . The minimum sampling 
rate for data to be used in this heading instrument (based on 
the design goal of Chapter II) then results in a system band- 
width of 30 hertz. The analog signals from each sensor are 
low pass filtered to reduce frequency content above 60 hertz. 

A survey of commercially available multiplexers, sample and 
hold modules and analog to digital convert modules (ADC) [Ref. 
3-25 to 3-28] reveals that subsystems with throughput -charac- 
teristics exceeding the requirements of a system sampled at 
one second intervals are readily available (pertinent specif- 
ications are discussed in more detail in Chapter IV) . The 
limiting parameter determining total system speed performance 
will ,then be the execution time of the algorithm (a programm- 
ing consideration) . A further system consideration is the 
ability to adjust analog system of-fset and gain. These ad- 
justments are made using variable resistors (trim pots) con- 
nected to appropriate leads on the sample and hold and analog 
to digital converter modules. 

2) The Central Processing Unit (CPU) 

The central processing unit (Fig. 3-2) is composed of the 
microprocessor (Signetics 2650) supported by peripheral logic 
elements (Fig. 3-2) . Design of this subsystem involved medium 


^That is, the sampled data derived from a sine wave of freq- 
uency f sampled at a frequency less than 2f can be fitted 
with sine waves of a frequency other than f. 
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and small scale integrated circuits using well-known [Ref, 

3-29 through 3-31] design techniques. To facilitate system 
development several features were included in the design of 
the GPU subsystem (features that would not necessarily be re- 
quired in a production instrument). These include: 

a) System reset, single step and normal run mode oper- 
ation controlled by switches and logic elements. 

b) An RS-232 teletype interface is included to enable 
manual intervention and development capability 
during program development. The program was dev- 
eloped by loading and executing instructions into 
the random access memory (RAM) under control of the 
PIPBUG^ program. 

3) The Memory Subsystem 

The memory subsystem (Fig. 3-3) was organized onto cards 
each with two thousand byte capability. In this manner system 
memory could easily be expanded (or reduced) in increments of ' 
2K bytes. The memory chips selected were organized as 256 
four bit words and feature pin for pin compatibility with 
commercially available random access (RAM) and programmable 
read only memory (PROM) chips. Program segments could then be 
developed in RAM and finally "burned" into PROM chips for a 
permanent, nonvolatile operation. In this manner the system 
development begins with IK bytes of memory devoted to the res- 
ident PIPBUG program (in ROM chips) with the remainder of mem- 
ory allocated as RAM for both program and scratch pad usage. 


^Signetics tradename for the 2650 resident, loader and monitor 
program. 
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Fig. 3-3. MEMORY SUBSYSTEM (1 Card) 























37 


As the program is developed, additional memory is added in 
increments of 2K bytes per card or 256 bytes on the card. 
Modifications to the program can be easily made using the 
PIPBUG program and teletype. 

4} The Output Subsystem 

For laboratory development the output subsystem of Fig. 
3-4 was designed to provide seven segment visual output of 
the aircraft heading with three significant digits displayed. 
To expediate the design cycle and to enhance system throughput 
rate, the outputs were designed as ports with latches and de- 
coder driver functions provided by hardware . In other applica 
tions a hardware/software tradeoff could be made with the data 
decoding and driving implemented using table lookup and multi- 
plexing controlled by the CPU. 

3-3 SOFTWARE DESIGN CONSIDERATIONS 

'The general purpose processor selected to implement the 
CPU was designed to implement programmed logic and to perform 
conventional computer operations. This heading instrument 
takes advantage of both areas. Since the instrument is act- 
ually a special purpose computer under control of a stored 
program, the functional specialization resides in the program 
rather than the hardware logic. Modifications can be made re- 
latively easily, satisfying the flexibility design goal of 
Chapter II. 

Having decided on the tentative hardware structure de- 
scribed in Section 3-2 above, the program development leading 
to the final listing in Appendix B proceeded as follows: 

1) Structured flow charts were developed depicting the total 
system operation as an ordered sequence of operations . Each 
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INITIALIZE 


DO forever 


main 


power on reset of all registers 
and subsystems 


compute the aircraft heading 


Fig. 3-5a. 


SYSTEM PROGRAM 


SAMP 


3 




^ ROTY ^ 


hvec ^ 


WICH ^ 


^ OUTA ^ 


sample all analog channels 

Compute the horizontal Hx field 

compute the horizontal Hy field 

Compute the horizontal field 
vector 

Compute heading 
Output the data 


Fig. 3-5b. 


SUBROUTINE "MAIN 





DO 5 TIMES 


SELECT AN ANALOG 

CHANNEL 

SAMPLE AND HOLD 

DATA 

CONVERT ANALOG TO 

DIGITAL 

DO TWICE 


STORE 8 BITS 
OF DATA 


CONVERT UNIPOLAR BINARY DATA 
TO SIGN MAGNITUDE FORMAT 


CHANGE THE SIGN OP 
Hx AND Hy SENSOR DATA 


DO '3 TIMES 


OFST 


CORRECT FOR SENSOR OFFSET 


ORTH 


CORRECT FOR SENSOR 
ORTHOGONALITY ERROR 


RETURN 


Fig. 3-6a. SUBROUTINE "SAMP 
















PRELOAD VARIABLES REQUIRED 
TO CORRECT FOR OFFSET ERROR 

^ SAPP ^ 

PERFORM A DOUBLE PRECISION ADD OR SUBTRACT 
PUSTORE CORRECTED DATA TO TABLE 
RETURN 

Fig. 3-6b. SUBROUTINE "OFST" 


TEST 

FLAG 

ADD 

SUBTRACT 

DO DOUBLE PRECISION 


ADD ON SIGN 

DOUBLE PRECISION 

MAGNITUDE NUMBERS 

SUBTRACT 

RETURN 


Fig. 3-6c. SUBROUTINE "SADD" 
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input ports 1 and 2 and output port 1 of Fig. 3-1 are allocat- 
ed. Control information is then passed to the peripheral mod- 
ule by writing control words to output port 1. Analog to dig- 
ital converter status and the 12 bit data field are sampled by 
reading input ports 1 and 2. 

Sensor outputs were biased at +2.5 Volts with transfer char- 
acteristics as depicted in Fig. 3-7a [Ref. 3-34]. The ADC 
selected for this laboratory instrument had a binary output 
data format related to analog input as shown in Fig. 3.7b 
[Ref. 3-35]. The second function of the sampling subroutine 
"SAMP" was to convert data from a unipolar binary format to a 
sign magnitude format. Since the total transfer function from 
sensor input to ADC output (Fig. 3-7a and b) indicates an off- 
set of 2.5 Volts or 1/2 the ADC output range, the sign magni- 
tude format can be generated as shown in Fig. 3-8. 



Fig. 3-7a. 


SENSOR TRANSFER CHARACTERISTIC 



SIGN BIT 



Fig. 3-8 CONVERSION OF DATA 
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The third function of the "SAMP" subroutine was to reverse the 
sign of the Hx and Hy data (to correct a test fixture pro- 
blem) and to correct for sensor offsets. Although analog sub- 
system offsets are corrected by adjusting either the sample 
and hold module or the ADC, the independent sensors themselves 
have offsets^. Offset errors for the laboratory instrument 
were compensated by determining the offset correction term for 
each sensor (method described in detail in Chapter V) and then 
either adding or subtracting the term to the respective data 
during the sample subroutine. By characterizing the sensor 
errors^, actual datum could be improved further during this 
step. 


The final function of the "SAMP" subroutine was to cor- 
rect for sensor orthogonality error (subroutine "ORTH"). Al- 
though the sensors were physically aligned and specified to 
have orthogonality characteristic (Ref. 3-34] less than +1 
degree relative to the base coordinates, this nonorthogonality 
contributes appreciably to total system error (see error an- 
alysis in Chapter IV) . The physical misalignment of the sens- 
ors was determined experimentally (Chapter V) and determined 
to be mainly a misalignment of sensor x in the x-y plane as 
illustrated in Fig, 3-9, 

The actual data measured with the x axis sensor is then relat- 
ed to the true Hx and Hy values as 

Hx^ = Hx Cos e - Hy Sin e , 


^With zero stimulus applied the sensors have a finite nonzero 
output. This error in the fluxgate magnetometer is a function 
of temperature, voltage and magnetic remanence in the sensor 
magnetics [Ref. 3-36]. 

^Sensor characteristics relating the temperature and power sup- 
ply coefficients of offset error and nonlinearity can be de- 
rived empirically . 


46 



Fig. 3-9 X AXIS NONORTHOGONALITY 


Using small angle approximations, we can solve for the de 
sired true value of Hx 


Hx^ ^ Hx - Hy Sin £ 


(3-la) 


Hx == Hx^ + Hy Sin £ 


(3-lb) 


By measuring e (Chapter V) and storing the angle as a con- 
stant, the X axis data was then restored using equation 3-lb 

above in subroutine "ORTH" . 

2) Subroutines ROTX and ROTY 

These subroutines compute arithmetic values for Hxh and 
Hyh of equation 2- 11b using sign magnitude quantities and 
table lookup to determine solutions for the transcendental 
functions. Subroutines "SADD" and "SMPY" are nested and used 
to perform double precision add and multiply as required. 
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3) Subroutine HVEC 

Following computation of the horizontal X and Y axis mag- 
netic vector, the subroutine "MAIN" calls subroutine "HVEC" to 
compute the square of the horizontal vector. Vectors Hx and 
Hy are squared by calling subroutine "SQU" then added, yield- 
ing H (HORIZONTAL) 

4) Subroutine WICK 

To compute heading, equation 2-1 (or a similar form) must 
be solved using the horizontal magnetic field vector and either 
the X or y axis horizontal field component. Although the 
square root operation implied in equation 2-1 could be imple- 
mented using a numerical technique [Ref. 3-37, 3-38], the 
computation time is decreased by using a table lookup method. 
Subroutine "WICH" (Fig. 3-10) compares the absolute magnitude 
of the two horizontal field vectors Hx and Hy to determine the 
relative heading of the aircraft^ with respect to the north- 
south' and east-west axes (Fig. 3-11). 



Fig. 3-10. SUBROUTINE "WICH" 


^If |Hx|<|Hy|, then an equation similar in form to 2-la must 
be used. 
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N 



Fig. 3-U. MAGNITUDES OF Hx AND Hy RELATED 
AIRCRAFT HEADING 

5) Subroutines COSY and SINY (Fig. 3-12, 3-13) 

Depending on the relative absolute magnitudes of Hx and 
Hy, either "COSY" or "SINY" is called to compute aircraft 
heading. These subroutines invoke subroutine "DIVI" to form 
the quotient of the axis vector squared and the horizontal 
field vector squared (a double precision operation) . Sub- 
routine "ANGL" is then called to perform an associative table 
lookup operation using successive approximation and inter- 
polation to complete the inverse cos squared operation. The 
double precision binary quantity is then converted to three 
’ digit binary coded decimal format (BCD) prior to computation 
of aircraft heading (subroutine "HDG"), 

The subroutine "SINY" of Fig. 3-12 includes a subtraction 
of the computed angle from 90 degrees following conversion to 
BCD format. This operation ensures that the angle passed to 







Fig. 3-13. 


SUBROUTINE ”SINY 
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the calling subroutine upon exiting either "SINY" or "COSY" is 
an aircraft heading angle relating sensor x to the north-south 
axis . 

6} Subroutine HDG (Fig. 3-14) 

The function of this subroutine is to compute aircraft 
heading having established the angle between the x axis sen- 
sor and .the north-south geodetic axis. Determination of the 
actual heading is accomplished by comparing the signs of both 
the x and y axis horizontal vectors prior to computing heading 
(Fig. 3-15). It should be noted that all of the preceding 
computations leading to horizontal vector data were on sign 
magnitude quantities preserving the correct horizontal vector 
polarities'^ . 


3-5 CONCLUSIONS 

This chapter has outlined the practical aspects of design- 
ing an instrument to evaluate both the heading algorithras and 
solid state magnetic indicator proposed in previous chapters. 
The chapter outlined a design approach that can be used to 
implement a microprocessor based instrument. In particular, 
the need to consider the total system hardware requirements 
while simultaneously considering the programming requirements 
was identified. Design proceeded by outlining a system block 
diagram (Fig. 3-1) with major subsystems considered. The in- 
strument required a special purpose computer with an analog 
subsystem to sample and digitize five sensor signals. Timing 
and control of the analog subsystem plus digital processing 
of data was controlled by a microprocessor based central pro- 
cessing unit (CPU) . Memory for permanent storage of 


^It is possible at certain attitudes to require sign revers- 
als when computing horizontal vectors , 
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instructions and temporary storage of data was implemented 
using memory chips organized on cards with 2048 byte capacity. 
The particular memory chips selected feature pin compatibil- 
ity® with both read only and volatile random access versions. 
System inputs consisted of sensor signals from a three axis 
solid state fluxgate magnetometer plus two analog signals sim- 
ulating gyroscope outputs. System outputs consist of visual 
seven segment readout displaying computed heading. In addi- 
tion, an RS-232 teletype interface was provided to facilitate 
system development and experimentation. 

By identifying the total system in block diagram form at 
the very beginning, the role and requirements of each subsys- 
tem as well as the supporting software were identified. The 
design then evolved on a modular basis with each subsystem and 
its supporting program developed in parallel. In this manner 
pin assignments for input/output ports and critical timing re- 
quirements that involved both hardware and software considera- 
tion were handled efficiently. By outlining the program re- 
quirements in flow chart form (analagous to the block diagram 
of the hardware subsystem) , subroutines were identified facil- 
itating a modular program development. Where possible, sub- 
routines were shared in a nested manner avoiding replication of 
programming and waste of memory. 

Details of error analysis and calculation of overall sys- 
tem throughput rate were deferred to Chapter IV. It was point- 
ed out however, that errors induced by imprecision of data plus 
truncation and roundoff during processing of the algorithm were 
to be considered early in the design phase. These data were 
required to select the sensors and the analog to digital 


Memory integrated circuit (IC) devices of both types can be 
used in the same mechanical sockets with actual chip type 
being used transparent to the remainder of the system. 
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converter as well as to design the supportive software for 
the analog subsystem. In addition, the data precision require- 
ments were necessary prior to programming the algorithm^ . 

By incorporating a microprocessor as the main CPU element, 
considerable sophistication in both control and computing per- 
formance was achieved. The overall system was designed rel- 
atively quickly, provided a convenient laboratory instrument 
for evaluation of the proposed algorithms and featured inher- 
ent flexibility. 


^Some of the subroutines required double precision manipulations 
to maintain overall system accuracy. 
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CHAPTER IV 

HEADING INSTRUMENT ERROR ANALYSIS 
4-1 INTRODUCTION 


The heading instrument designed to evaluate the heading 
and solid state remote magnetic indicator algorithms is prone 
to error from many sources. These errors will accumulate and 
degrade the accuracy of aircraft heading or yaw angle computa- 
tions. This chapter addresses the various error sources to 
determine their relative magnitudes and effects on the overall 
computation. 

Prior to beginning the hardware design of the micropro- 
cessor based instrument many of these potential error sources 
were considered. Their effects were considered in, establish- 
ing parameters such as word lengths, A/D converter precision, 
computation speeds, sampling rates, magnetometer sensor accur- 
acies, system noise tolerance, etc. As the design of the 
microprocessor based system evolved, the error analysis re- 
fined. Ultimately, important limitations in instrument design 
and operation were identified by combined error analysis and 
empirical data. By carefully analyzing the source and extent 
of the limiting parameters (such as sensor offset and non- 
orthogonality) , the magnitude of errors unique to this lab- 
oratory sensor array were identified. Specialized software 
was then added (with empiracally derived constants) to correct 
for the otherwise limiting sensor irregularities improving the 
total system performance. 

In this manner, it is apparent that error analysis is an 
integral part of instrument design. Not only are important para- 
meters identified early in the design cycle (prior to system 
block diagram development) , but shortcomings in conventional 
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•sensors can be improved by judicial application of error cor- 
recting algorithms. In this case, data constants were deter- 
mined after the final instrument became operational. The sen- 
sor peculiarities were analyzed empirically using the instru- 
ment itself. 

The chapter begins by first identifying and carefully 
analyzing potential error sources in the sensors. This analy- 
sis is followed by a similar consideration of errors origin- 
ating in the analog subsystem. Processing errors that origin- 
ate due to the finite word length and precision of the micro- 
processor along with the effects of simplifications made to 
the algorithms are finally analyzed. The chapter then con- 
cludes with a summary of measurement errors, a sample error 
analysis, a comparison of predicted to measured error and a 
summary. 

4-2 SENSOR ERRORS 

The heading computation algorithm employing the remote 
magnetic indicator {Chapter II) is prone to error proportional 
to both fluxgate magnetometer sensor and gyroscope measurement 
errors. Errors inherent in the fluxgate magnetometer are sum- 
marized on the data sheet [Ref. 3-34]. Since the experiment- 
ation employed simulated gyroscope sensors with voltage levels 
accurately represented, the analysis of sensor errors will 
assume ideal gyroscope sensors to predict experimental data. 

A) Sensor Offset Error 

Magnetometer sensors exhibit error caused by both elec- 
tronic and magnetic phenomena. Errors in the Develco sensors 
were outlined by Workentine [Ref. 4-1] . These offset errors 
are induced in the Develco sensors by both electronic offset 
voltages and currents in the respective sensor electronics and 
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by residual magnetic fields in the magnetic mass of the sensor 
assemblies. Although the physical and electronic design at- 
tempts to reduce offset error, a finite non-zero output can 
exist when a zero input is applied. 

Offset error for each sensor in the Develco model 9200C 
three axis magnetometer assembly is specified [Ref. 3-341 as 
"Zero Field Bias 4-2.5 Volts +1.0%**. This offset translates 
into a worst case maxim\im error voltage of 

EOFFSET = +{2.5V X 0.01) = +25mV 

Since the offset error is sensor dependent, correction cannot 
be made at a single physical point (as for analog subsystem 
offsets described in Section 4-3) . Corrections can however be 
made to the measured data by simply adding or subtracting a 
constant equal to the offset magnitude following each data 
measurement^ . 

Offset values for each sensor used in the experiment were 
obtained by rotating the sensor into alignment with earth's 
magnetic field vector to measure both positive and negative 
maximum values. The difference in magnetic measurement (assum- 
ing negligible analog subsystem error) is related to system 
offset error composed of sensor electronic and sensor plus test 
fixture induced magnetic offset error. The actual offset error 
can be calculated using these two measurements 


^Offset corrections were made in the sample subroutine "SAMP 
illustrated in Fig. 3-6a, 
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|Emax| = Ef + Eo 
I Emin I = Ef - Eo 

Emax| - I Emin I = (Ef + Eo) - (Ef - Eo) = 2 Eo 
Eo = (1/2) (|Emax) - | Emin | ) 


where 

Emax The maximum positive voltage recorded when the 
sensor aligns with earth’s field vector. 

Emin The maximum negative voltage recorded when the 
sensor aligns 180* with earth's field vector. 

Ef = The magnitude of earth's magnetic vector repre- 
sented in volts. 

Eo = The sensor offset voltage due to both electronic 
and magnetic phenomena 

Data recorded during x, y and z axis offset measurements as 
described above are recorded in Table 4-1. Since the offset 
error is a function of sensor magnetic permeability, the actual 
offset value will vary with time depending on induced magnetic 
fields^ . 

Final offset correction values were determined by rotat- 
ing two sensors in the horizontal plane around the third ver- 
tical axis and measuring offsets in two sensors at a time. 
Recorded data for each sensor was previously corrected for 
orthogonality error by the sample subroutine "SAMP" (discussion 


2por example, magnetized screwdrivers or other tools used near 
the sensor will alter the residual magnetic field. 
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SENSOR 

AXIS 

DATA RECORDED 
1 Eroax 1 

(HEXADECIMAL) 
1 Emin 1 

OFFSET 
(HEXADECIMAL 
PLUS SIGN) 

X 

628 

E78 

+ 40 

Y , 

640 

E5D 

-15 

Z 

637 

E68 

-25 ■ 


Table 4-1 OFFSET DATA DERIVED BY MEASURING 
EARTH'S FIELD 
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of this correction follows in Section 4-2B) . Data recorded in 
this manner appears in Tables 4-2 and 4-3. Final correction . 
terms for correcting sensor offset error were calculated using 
these data. Offset terms to be added or subtracted from res- 
pective data channels are tabulated in Table 4-4. 

By correcting system offset errors in this manner, the 
effective error contribution can be reduced appreciably (see 
final data discussion Chapter V) . For a flight instrument, 
sensor offset characteristics as a function of temperature 
variation and supply voltage can be derived empirically and 
appropriate offset corrections made by computing the value of 
the correction term variable. Magnetically induced offsets can 
be reduced by degaussing the sensor assembly periodically. 

B) Axis Alignment Errors 

The error specification of [Ref. 3-34] indicates that the 
maximum axis alignment error is +1 degree relative to base ref- 
erenced coordinates. This error results in sensor directional 
uncertainty as illustrated in Fig. 4-1. Each sensor is located 
within a right circular cone with axis along the true sensor 
axis and vertex at the common sensor origin. Although this 
alignment uncertainty contributes no error in determining the 
total magnetic vector 


_ _2 _2 _2i 

H = (Hx + Hy + Hz )^, 

there is considerable uncertainty in attempting to resolve the 
true magnetic field component along any axis of the reference 
coordinate system. This alignment uncertainty of magnetic 
sensors limits system performance of conventional field dir- 
ection measuring apparatus [Ref. 4-1]. 
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protractor 

Heading 

Measurement 

(Degrees) 

0 

345 
330 
315 
300 
285 
270 
255 
240 
225 
210 
195 
180 
165 


Data 

Measured 

Hx Hy 
(Units) 

8 -759 
199 -743 
383 -677 
540 -569 
666 “421 
749 -249 
783 - 62 
770 132 

706 314 

597 473 

449 601 

273 687 


Protractor Data 

Heading Measured 

Measurement Hx Hy 

(Degrees) (Unrts) 


180 

165 

150 

135 

120 

105 


83 720 

■112 704 

-292 638 

-448 527 

-573 381 

-658 206 


Error Due To 
Offset 
X Y 

(units) 


75 

87 

91 

92 

93 
91 


90 

-696 

14 

87 

75 

-681 

-179 

89 

60 

-616 

-362 

90 

45 

-504 

-521 

93 

30 

-359 

-645 

90 




■39 

-39 

-39 

-42 

-40 

-43 

-48 

-47 

-48 

-48 

-44 


total offsets 


978 


477 


average offsets 


88.9 43.4 


Table 4-2 


Mn Y AXIS ERROR MEASURED BY ROTATING 
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Protractor 

Heading 

Measurement 

(Degrees) 

Hz Data 
Measured 
(Units) 

Protractor 

Heading 

Measurement 

(Degrees) 

Hz Data 
Measured 
(Units) 

Offset 

Error 

(Units) 


0.5 

Q 

180.5 

-56 

-56 


315.5 

-527 

. 135.5 

469 

-58 


270.5 

-763 

90.5 

707 

-56 


225.5 

-571 

45.5 

513 

-58 



TOTAL OFFSET 


-228 



AVERAGE OFFSET 


-57 


Table 4-3 

Z AXIS 

OFFSET ERROR MEASURED BY ROTATING 



THE Z AXIS AROUND THE 

VERTICAL 

X AXIS 


Sensor Total 

. Average 

Required 

Amount ( 

Df Correction 

Axis Offset (Units) 

Correction 

Decimal 

Binary 

Hex 

X 

88.9 

Subtraction 

45 

00101101 

02D0 

V 

43.4 

Addition 

22 

00010110 

0160 

s 

47.0 

Addition 

29 

OOOUlOl 

OIDO 


Table 4-4 OFFSET CORRECTION VALUES 



62 



Fig. 4-1 SENSOR ALIGNMENT UNCERTAINTY 

Although this error source can be reduced by physically 
aligning the sensors more accurately during assembly, cost of 
the sensors increases. Ultimately, directionality of the mag- 
netic sensors becomes a function of the physical sensor itself 
and more accurate sensors are required as pointed out by Gise 
[Ref 4-2]. A heading system that tolerates sensor misalign- 
ment' is therefore a very desirable alternative to requiring 
precise alignment or more elaborate sensors. 

During assembly of the Develco fluxgate magnetometer sen- 
sor array, sensor misalignment is determined by using earth’s 
magnetic field and a precision mechanical rotation assembly. 

A sensor (assume the X axis) is aligned with earth's magnetic 
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vector by positioning the sensor to maximize electrical output^ . 
One of the other sensors (assume the y axis) is aligned with the 
rotation axis of the precision calibration assembly (Fig, 4-2) 
and perpendicular to the first by rotating the sensor array 
around the second sensor axis (y axis in this case) and adjust-, 
ing its relative position until a null output is achieved at all 
rotation angles. Mechanical orthogonality of the sensors is then 
limited only by the mechanical imprecision of the calibration 
device (orthogonality within +0.01 degrees can be easily achieved 
in the calibration tool) and by the directional characteristics 
of the physical sensors. 


X Axis 



Fig. 4-2 MECHANICAL ORIENTATION OF THE MAGNETOMETER 
SENSORS DURING CALIBRATION 


In addition to functioning as an alignment apparatus, the 
calibration device described above provides a convenient means 


^By maximizing or nulling a measurement, the mechanical posi- 
tioning is a function of only the field and the resolution of 
the voltage measuring device obviating errors due to physical 
position measurement. 
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to characterize sensor assemblies after final assembly adjust- 
ments are made. Any misalignment of the second sensor relat- 
ive to the first results in a coning of the second sensor 
around the rotation axis'* with a sinusoidal output voltage 
that is a function of total earth's magnetic field and axis 
alignment error. The peak to peak voltage resulting from sen- 
sor coning is recorded during the final alignment test and 
made available to -sensor purchasers. Coning voltages developed 
for the sensor assembly used with this experiment were obtain- 
ed from Develco [Ref, 4-3] and are recorded in Table 4-5, Sen- 
sor misalignment for each axis can be derived using additional 
data provided by Develco along with additional empirical data 
derived by experimentation. 

The total ambient magnetic field at the Develco laboratory 
is measured using the three sensors (applying equation 4-1) 
and is supplied as digital data. In our case, the total field 
measured was 1573 units or 


1573 units 

2048 units 
Full Scale 


X 60,000 gamma F.S. = 46,084 gamma (y) 
(F.S.) 


Sensitivity of the sensor = - ^2 uVolts/y 

oU , {jUU Yr . b . 

Considering the X axis sensor, coning resulted in a signal of 
38 mV peak to peak (or 19mV peak) . Misalignment of the X axis 
sensor from the Y-Z plane can then be calculated as 


‘^Assume that the first axis is initially adjusted for maximum 
output to align it with earth's field and the rotation axis 
is perpendicular to the field. 
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Rotation 

Axis 


X 

Y 

Z 


SENSOR ASSEMBLY NO. S/N 1043-013 


Coning Voltage Orthogonality 

(Peak-Peak mV) Error (Degrees) 


38 


0.57 


8 


=*0 


51 


0.76 


Table 4-5 MAGNETOMETER ORTHOGONALITY MEASUREMENTS 
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Peak Signal = 19 mV or 456 gamma angular misalignment 

= Sin ^ 456 

^ 46,084 

= 0.57 degrees 

Similarly, the Y and Z axis have misalignment errors of 
ey 0 and tz = 0.76 degrees with respect to the X-Z and X-Y 
planes respectively (sensor orthogonality errors are tabulated 
in Table 4-5} . 

Having established that sensor orthogonality errors exist, 
the remaining task is to identify the direction that the sensor 
axis points relative to the other two sensor axes. Since the 

Y axis has relatively little orthogonality error, it will be 
assumed to be perpendicular to the X-Z plane. In addition, 
since the Hz data enters into the algorithm in a second order 
manner relative to the Hx and Hy measured data, correction and 
characterization of the Hx sensor was considered to be of prim- 
ary concern. Orientation of the X axis sensor relative to the 

Y and Z axes was determined empirically. 

Angular position of the X axis sensor can be described 
using the error angles exy and exz as delineated in Fig. 4-3. 
Characterization of sensor orthogonality error in terms of 
these two angles would enable algorithmic corrections of mea- 
sured data. 



Fig. 4-3 X AXIS SENSOR ORIENTATION 
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1) Empirical Determination of txz 


::;t ::.:rr;r. ’■ 

X .nd , ..«or, o.l»t.a in «.• h«U=«l.l 

pL. dim tn. , .ni. 


i) 


ii' 


iii) 


..e X an. . seniors .ere L^reL.. 

with magnetic data measurements (co 

offset error as described rn Sectron 4 2A r 
in Table 4-6 for incremental rotation ang e • 

..e total horizontal field at each angular position 
was calculated 

Hht = (Hx2 + Hy2)^ 

f-sl field Hav was computed by averag- 
Average horizontal fieia nav 

ing the results of iii) above. 

U 'anneal field deviation Hd was computed for 
The horizontal fie Table 4-6 and 

each angular position; tabulated in 

plotted on Fig- 4-4. 

Hd = (Hav - Hht) 

n„nonn.l ll.ld deviation .n .»on 

ananlned. » d.,i=a. ^ 

horizontal field error curv concluded that 

Since this obviously was not 

error in x axis orthogonality was 

£xy . 


i'v) 


V 
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Physical* Displayed** 
Heading Heading 

(Degrees) (Degrees) 


Total 

Computed 

Measured Data Horizontal 
(Units) Field (Hht) 

(Units) 


355 
335 
315 
295 
275 
255 
235 
215 
195 
175 
155 
135 
115 
95 
' 75 
55 
35 


90 

12 

70 

258 

50 

477 

30 

637 

10 

720 

350 

718 

330 

628 

310 

464 

290 

242 

270 

-10 

250 

-258 

230 

-476 

210 

-635 

190 

-721 

170 

-718 

150 

-628 

130 

-460 


-727 

727 

-682 

729 

-555 

732 

-362 

733 

-127 

731 

127 

729 

363 

725 

569 

734 

687 

728 

731 

731 

686 

733 

557 

733 

362 

731 

123 

731 

-133 

730 

-371 

729 

-565 

729 


Hd 

(Hav-Hht) 

(Units) 



Total Hht = 12415 


Average (Hav) = 730 


^Measured using a protractor on the test 
**Computed and displayed digitally by the 


apparatus . 
instrument . 


Table 4-6 MEASUREMENT OF HORIZONTAL FIELD 



Error (Units) 



Fig. 4-4 deviations 
from the M 



Error (Units 



Fig. 4-5 DEVIATION OF Hx AND Hy • DATA FROM THE COMPUTED 
FIELD COMPONENTS AS A FUNCTION OF YAW 
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2) Empirical Determination of exy 

The angle exy representing x axis sensor misalignment 
relative to axis y was measured as follows: 

Steps i) and ii) above were repeated with the exception 
that the calculated values for Hx and Hy (Hxc and Hyc respect- 
ively) were recorded with measured Hx and Hy data (Hxm and Hym 
respectively) in Table 4-7. The calculated values were obtain- 
ed by assuming that the angle exz as determined above was neg- 
ligible and that the y axis sensor was perpendicular to the x-z 
plane. With these assumptions, we note that at the heading of 
zero degrees (extrapolated between display of 10 and 350 degrees 
of Table 4-7 and Fig. 4-5), there is no error in yaw due to 
either Hx or Hy. By physically rotating the sensors in fixed 
intervals from yaw = 0 degrees and noting that the horizontal 
field Hh = 730 units, we can then compute expected Hx and Hy 
data at respective yaw orientations. 

Physical orientation of the x axis sensor is easily deter- 
mined by considering orientation at the maximum epror excurs- 
ions. These observations are illustrated in Fig. 4-6. We note 
that the only possible orientation of the x axis sensor satis- 
fying the data in Fig. 4-5 is that of Fig, 4-6. 



Fig. 4-6 (a) SENSORS ORIENTED AT YAW = +90 degrees 
(b) SENSORS ORIENTED AT YAW = +270 degrees 


Physical* Displayed** Measured Data Computed Data 
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Table 4-7 MEASURED AND COMPUTED Hx AND Hy DATA 
IN THE HORIZONTAL PLANE 



73 


Magnitude of the angle exy can be computed as follows 
using data from Fig. 4-5 

Max. delta from Fig. 4-5 == 10 units 
Average horizontal field =5 730 units 

exy max = Sin ^ 

= 0.79 degrees 

We note that the angle of 0.79 degrees is approximately 
the same as determined by Develco during manufacture of the 
sensors (Table 4-5) . The added error is due to test set in- 
accuracy. 

C) Fluxgate Sensor Noise Induced Error 

The analog output from the fluxgate sensors can exhibit an 
error due to signal uncertainty resulting from noise. Although 
the data sheet [Ref. 3-34] indicates that 5mV peak to peak of 
ripple can exist on the output, the frequency content centers 
in the 550 kHz range (driver frequency of the fluxgate magneto- 
meter) and no appreciable ripple^ exists below 60 Hz (especial- 
ly when the sensor output is filtered prior to data sampling) . 
The noise specification of less than 1 gamma peak to peak in 
the 1' Hz bandwidth region is also negligible. In summary, no 
appreciable error due to noise on the magnetometer signal lines 
is evident. 

D) Magnetometer Gain Error 

The magnetometer is specified to have gain (sensitivity) 
of 2.5 Volts/600 milligauss, +1% which translates into a maxi- 
mum signal uncertainty of 

^Verbally confirmed by Workentine of Develco [Ref. 4-1]. 
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+(2.5V X 0.01) = +25 mV. 

Chis represents a sensor transfer function of 4.16 Volts/gauss 
Dr 0.24 gauss per volt. The uncertainty then can be expressed 
as 


+(0.24 gauss x 0.01) = +2.4 milligauss 

= +(2.4 X 10^) gamma 

Since this error is not corrected in the laboratory instrument 
it will be considered in total in the final error analysis. It 
is worth noting- however, that should the magnetometer gain un- 
certainty be characterized, gain corrections for each sensor 
could be made during computation by the computer. In addition 
the error term is proportional to actual signal level applied. 

E) Magnetometer Linearity Error 

D.C. linearity of the magnetometer is specified to be 
+0.5% of signal level. This uncertainty at full scale can be 
expressed as +(2.5 Volts x 0,005) = +l2.5mV, Alternately, lin- 
earity error can cause a signal uncertainty of +1.2 milligauss 
or +(1.2 X 10^) gamma. Linearity error is also not corrected 
during computation and is considered in the final error analysi 
By simply characterizing and correcting the linearity character 
istics of each sensor, considerable improvement in system accur 
acy could be achieved. 

4-3 ANALOG SUBSYSTEM ERROR ANALYSIS 

•The analog subsystem of the instrument is outlined in 
block diagram form in Fig. 4-7. This subsystem accepts analog 
signals from magnetometer and gyroscope transducers , performs a 
time division multiplexing between the signals and digitizes 



-15 +5 +15 



Fig. 4-7 THE AtlALOG SUBSYSTElyi 
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the respective signals prior to subsequent processing, by the 
computer. During this data acquisition and conversion process, 
errors are introduced into each of the signals. This section 
addresses the potential error sources and computes the respect- 
ive error contributions to be expected during operation of the 
instrument. 

Although the multiplexer and sample and hold blocks of Fig. 
4-7 could be eliminated (eliminating possible error sources) by 
digitizing each signal with a unique analog to digital converter, 
it can be shown that such a system would be expensive and dif- 
ficult to implement. The analog to' digital converter (A/D) 
quantizes an analog signal in a finite amount of time. Speed of 
conversion is predicted in a finite amount of time by both the 
resolution of the converter and the frequency of the signal to 
be converted. Time required to perform a conversion is general- 
ly called 'the "aperature time". 



Fig, 4-8 APERATURE TIME AND AMPLITUDE UNCERTAINTY 

As illustrated in Fig. 4-8, aperature time and amplitude 
uncertainty are related by the time rate of change of the ana- 
log signal. For the particular case of a sinusoidal signal to 
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be converted, the maximum rate of change occurs at the zero cross- 
ing of the waveform and the amplitude change is : 


AV = ^ (V Sin wt)^ ^ Q X ta (4-la) 

AV = V w ta (4-lb) 

giving ^ = w ta = 2 irf ta. (4-2) 

From this result we can determine the aperature time required to 
digitize a 30 Hz signal to 12 bits resolution (a resolution of 1 
part in 4096 or 0.0244%). 


_ V 1 _ .000244 

V ^ 2 f 6.28 X 30 


1.3 X 10"® 


This result indicates that to remain within 1 bit of resolution 
(0.0244%) we require an aperature time of 1.3 microseconds to pro- 
cess analog signals varying at a rate of 30 Hertz. It can be seen 
that the system would require fast A/D converters plus extremely 
fast computational capability to accommodate this configuration of 
sensors and analog subsystem. By using multiplexing and sample 
and hold circuitry we can however reduce the number of A/D con- 
verters required to one and alleviate the aperature and processing 
requirements imposed above. 

The operation of sampling to be used by the instrument is 
illustrated in Fig. 4-9 which shows an analog signal and .a train 
of sampling pulses. The pulses are provided by the central pro- 
cessing unit. A switch connects the analog signal for a very 
short period of time to the hold circuitry charging a capacitor 
and storing the sampled voltage until the next sample is required. 
This type of sampler is called sample and hold. 
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A) Sampling Rate Errors 

The process of uniformly sampling a function of contin- 
uous time can yield a significant source of error if the sampl- 
ing period T is selected too large [Ref. 4-4, 4-5]. ' This error 
can be illustrated by considering an analog signal xa(t) that 
has the Fourier representation [Ref. 4-6] 


xa(t) = y__^Xa C jfi) e^^^dQ (4-3a) 

Xa{jQ) - r *xa(t)e ^^^dt (4-3b) 


The sequence x(n) with values x(n) = xa(nT) is said to be de- 
rived from xa(t) by periodic sampling and T is the sampling 
period. The reciprocal of T is called the sampling frequency 
or sampling rate. In order to determine the sense in which 
x(n) represents the original signal xa(t), it is convenient to 
relate Xa(jfi), the continuous-time Fourier transform of xa(t); 
to X(e^^), the discrete-time Fourier transform of the sequence 
x{n). From (4-3a) we note that 

x{n) = xa(nt) = ~ /'_"xa{ jn) (4-4) 


From the discrete-time Fourier transform we also obtain the 
representation [Ref. 4-4] 

x(n) = ^ r ’'X{e^“)e^“"dai (4-5) 

2k J-i> 


To relate the equations (4-4) and (4-5) we can express (4-4) as 
a sum of integrals over intervals of length 2 tt/T, as in 



{2r+l)7r/T . 

Xa(jn)e^“ dfi 


(4-6) 
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LCh term in the sum can be reduced to an integral over the 
inge -tt/T to +tt/T by a change of variables to obtain 


to 

r tt/t 

[- 

X (n) = — ^ XI 

/ Xa 

J -7t/T 

-i ( 

x(nj 2tt ^ 

r=-« 







(4-7a) 


X (n) 



^ .2-nr. ji^nT j2Trrn, 

Xa(3Q+3-^)e-' e-* dQ 


(4-7b) 


: we now change the order of integration and summation and 
)te that = i for all integer values of r and n, we ob- 

lin 
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substituting ^ = w/T we get 


X n = 


2 ir 



CO 

/ 1 

iE xa(^j^) 

J -TT. 

r=-o3 




(4-9) 


lich is identical in form to equation (4-5) . We can therefore 
ike the identification (equating like terras of (4-5) and (4-9) 


X(e^“) = i E Xa(^j^) 
r— -oo 


(4-10) 


i can also express (4-10) in terms of the analog frequency var- 
ible Q (where Q = oj/T) as 


X(e^“'^) = ^ E Xa(jfi+j^) 
r=-“» 


(4-11) 
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The last two equations clearly reveal the relationship 
between the continuous-time Fourier transform and the Fourier 
transform of a sequence derived by sampling. For example, if 
Xa (jn) is as depicted in Fig. 4-lOa then X(e^'^) will be as 
shown in Fig. 4-lOb when the sampling period T is too long and 
as shown in Fig. 4-lOC if T is short enough. 

QoT 

From Fig. 4-lOc it is obvious that if - <it, i.e., we 
sample at a rate at least twice the highest frequency of Xa(jfl) , 
then X(e^^) is identical to Xa(w/T) in the interval and 

can be recovered from the samples xa(nT) by an appropriate 
interpolation formula . 

For the remote magnetic indicator instrument designed in 
previous chapters, the analog signals are filtered with a low 
pass section reducing frequency content above 30 -Hz. The sampl- 
ing rate must therefore exceed 60 Hz (T<16.67 m.s..) to enable 
accurate dynamic operation of the system. 

Laboratory measurements of sampling rates on the function- 
al microprocessor based instrument revealed that the analog 
subsystem operated at a sampling rate of 62.5 Hz (16 m.s.) in- 
dicating that the algorithm execution rate supported a system 
bandwidth of 31.25 Hz. If frequency content of the analog 
signals is less than 31.25 Hz there is no error due to sampl- 
ing . 

B) Analog Multiplexer Induced Error 

The analog multiplexer of Fig. 4-7 selectively connects 
one analog transducer output at a time to the input of the 
sample and hold subsystem. The Datel Systems, Inc., multiplexer 
[Ref. 4-7] selected for the remote magnetic indicator experiment 
features eight MOS-FET switches with associated driver circuits. 
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FET pull-up to reduce propogation delays and all of the nec- 
essary decoding logic to enable random channel addressing with 
a four bit parallel binary input. 

Several important parameters are used to characterize ana- 
log multiplexers and can contribute error. 

1) Transfer Accuracy 

Transfer accuracy is a function of the source impedance, 
switch resistance, load impedance {if the multiplexer is not 
buffered) and the signal frequency. It expresses the input to 
output error as a percentage of the input. In our case the 
system configuration predicates a maximum error due to transfer 
accuracy of (+0.01%) yielding an error term of 

+0.0001 X 2.5 Volts ^ +25 mV 


2) Settling Time 

This parameter defines the time elapsed from the applica- 
tion of a full scale step input to the time when the output has 
entered and remained within a specified error band around its 
final value. In our case the selected multiplexer has a maxi- 
mum settling time of 1 microsecond to +0.01% full scale (F.S.) 
Since the control system selecting channels is implemented 
using a microprocessor, the minimum time between analog sub- 
system commands will always be greater than 3.0 'microsecond^. 
The multiplexer will therefore always have settled to the final 
value before the sample and hold circuit (following this sub- 
system) can be activated with no error due to the settling time 
parameter. 


^One machine cycle, time for the 2650 microprocessor with 1 mHz 
clock frequency. 
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3) Throughput Rate 

The highest rate at which the multiplexer can switch from 
channel to channel at its specified accuracy is in this case 
500 kHz. Since this rate is more than four orders of magnitude 
greater than the operational rate of the subsystem there is no 
error due to throughput rate limitations. 

4) Input Leakage Current 

The amount of signal coupled to the output as a percentage 
of input signal applied to all OFF channels together can be cal- 
culated by considering the maximum leakage current specified 
from OFF channels to the ON channel. In our case the maximum 
error signal can be calculated 

Error = [ 4 (8 na x 2000 ohms source imped. 

Error = 32 microvolts 

Note that in this case the voltage levels are statistically 
independent allowing an R.S.S. of error sources to calculate 
total error [Ref. 4-8, 4-9]. 

C) Sample and Hold Circuit Induced Errors 

The sample and hold subsystem consists of a switch and 
capacitor arrangement as shown in Fig, 4-11. The Datel Systems, 
Inc., model SHM-IC-1 integrated circuit sample and hold device 
[Ref. 4-10] features a self '■contained high gain differential 
input amplifier, a digitally controlled electronic switch and 
a high input impedance buffer amplifier. The external compon- 
ents used with the sample and hold circuit in the solid state 
remote -magnetic indicator instrument consisted of the O.OOlpf 
holding capacitor and a lOOK offset trimpot. By connecting 
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the output back to the negative input of the input amplifier 
(Fig. 4-11), the sample and hold subsystem operated in a unity 
gain, noninverting mode. When the switch is closed, the unit 
is in the sampling or tracking mode (Digital Control = 0 Volts 
and will follow a changing input signal. 



Figure 4-11 SAMPLE AND HOLD SUBSYSTEM 
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When the switch opens the unit is in the hold mode and 
retains a voltage on the capacitor for some period of time 
depending on capacitor and switch leakage. - Sample and hold 
devices are characterized by a number of important parameters 
that must be considered in the design of a data acquisition 
subsystem. 

1) Acquisition Time 

The time lapse between the time that the sample command 
is given to the point where the output enters and remains with- 
in a specified error band around the input value is specified 
to be less than 4 microseconds time to transit from 0 to 0.1% 
of 10 Volts with C = 0.001 jjf [Ref. 4-10]. This implies that 
the control signals emanating from the central processor should 
allow at least 4 ps acquisition time prior to entering the hold 
mode. We note that the sample and hold subroutine {Appendix 
B) executes the instruction 

lORI, R3 H'BO' READY TO HOLD, 

a two machine cycle instruction prior to sending the hold con- 
trol signal. This instruction delays control signal transmis- 
sion by (2x3 us) = 6 us allowing the sample and hold circuit 
ample time to settle with no appreciable error due to the 
acquisition time parameter. 

2) Hold Mode Voltage Droop 

The maximum change in output voltage as a function of time 
is specified to be 50 mv/sec maximum using a 0.001 uf poly- 
styrene capacitor. Since the maximum total accumulated time to 
completion of the analog to digital conversion can be calculated 


as 
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5 Instructions (11 machine cycles) 33 us 
1 Analog to Digital Conversion = 20 \is 

3 Instructions if Conversion not 
synchronized with instructions 
(7 machine cycles) = ^ ys 

74 ys 

we can then compute droop error to be 50 mv/sec x (74 x 10 ) 

sec — 3.73 mv . 

3) Aperature Delay 


The maximum time lapse between the time of hold signal 
receipt to opening of the switch is specified to be 50 nsec, 
an insignificant length of time in the instrument. There is 
therefore no error due to aperature delay. 

4) Offset Error 

Although the maximum offset error is specified to be 20 mv 
maximum (Ref. 4-10], the error was eliminated using the- lOOK 
trimpot offset adjustment. There was no appreciable offset 
error contribution due to the sample and hold circuit. 

5) Gain Error 

The gain error of a sample and hold circuit is apparent 
during the sample mode when the transfer function of the total 
amplifier deviates from the ideal unity slope condition (Fig. 
4-12) . In the noninverting unity gain mode, the specified gain 
error is +0.05% maximum yielding a signal error of 

(+0.0005) X 5.0 V = 250 mV 

This error can,, however, be eliminated with the gain 
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Fig, 4-12 GAIN, OFFSET AND LINEARITY ERRORS 

adjustment available at the analog to digital converter. There 
will 'therefore be no appreciable net gain error due to the an- 
alog subsystem. 

6) Nonlinearity Error 

Nonlinearity error is apparent in the sample and hold 
circuit if the transfer function departs from a linear curve 
(Fig. 4-12) . In the noninverting unity gain mode with a 
0.001 iif holding capacitor the maximum nonlinearity is 0.01% 
resulting in a worst case signal uncertainty of (.0001) x 
2.5 Volts = 25 mV. 

7) Hold Mode Feedthrough 

This error appears due to input signal appearing at the 
output when the unit is in the hold mode. Although the 
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feedthrough varies with signal frequency and the expected signal 
frequencies are substantially lower than the upper frequency 
limits of the sample and hold device (30 Hz max. versus several 
kiloHertz) , we consider the worst case feedthrough of 0.01% 

[Ref. 4-10] or 25 mV. . . 

D) Analog to Digital Converter Induced Errors 

The A/D Converter selected for the solid state magnetic 
indicator instrument (Datel ADC-MA12B1B) [Ref. 4-11] uses the 
successive approximation technique to achieve excellent lin- 
earity and speed. Important parameters that potentially con- 
tribute errors are addressed below. 

1) Resolution Error 


The smallest analog change that can be distinguished by 
the A/D converter is 


Least Significant Bit (LSB) 


Full Scale 
2 ^^ 


LSB = = 1.22 mV 

212 


this uncertainty manifests itself as an error in computing by 
limiting the precision of any calculation. 

2) Linearity Error 

The maximum deviation from a straight line drawn between 
the end points of the converter transfer function are specified 
in [Ref. 4-11] to be + 1/2 LSB (in our case +1,22 mV of analog 
signal) , 
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3) Accuracy Error 

The input to output error of the A/D converter is specif- 
ied in [Ref. 4-11] . to be +0.-012% F.S. +1/2 LSB or 

+{0.00012) X 5.00V +1.22 mV = +1.82 mV Worst Case 

In reality, the two error terms are unrelated and the 

Rss Error = + [(.00012 x 5V) ^ + (1.22mV)2]^ 

RSS Error = +1.36 mV 

4) Offset Error and Gain Error 

Both the offset error and gain error were adjusted to zero 
using the trimming potentiometers (Fig. 4-7) and the calibra- 
tion procedure outlined in Ref. 4-11. A reference signal of 
plus 1/2 LSB (1.22 mV) was applied to the converter and the 
offset trimming potentiometer adjusted until the output flick- 
ered equally between logic ”0" and logic "I". The gain was 
then adjusted by setting the converter input to full scale 
minus 1-1/2 LSB (4.99817 Volts) and the gain trimming potentio- 
meter was adjusted until the output flickered between logic 
"111... 110" and logic "111... 111". The above steps were re- 
peated until no appreciable error in gain or offset was evi- 
dent. 


4-4 PROCESSING ERRORS 

Errors in processing data accrue due to several sources 
including imprecision and truncation. Since the microprocessor 
selected for the instrument is inherently an eight bit device, 
single precision calculations are conducted with eight bits 
and double precision calculations are conducted with a total 
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of sixteen bits. This section addresses the effects of com- 
putational precision and truncation in the various subroutines 
and relates these to overall computational accuracy. The var- 
ious subroutines are analyzed in chronological order as they 
appear in the main program. 

A) Subroutine "SAMP" 

The sample subroutine (delineated in Fig. 3-6a) selects 
and digitizes analog signals by controlling respective analog 
subsystem modules. During the first portion of this subrout- 
ine, A/D converter data bits are stored in two consecutive 
bytes^ in the computer memory. The A/D conversion precision 
of 12 bits is thereby preserved. 

The second, third and fourth operations of the sample sub- 
routine convert the unipolar binary format of the data to sign 
magnitude format, adds offset quantities and merely changes the 
signs of the Hx and Hy data. The operations are conducted in a 
double precision manner and precision of the data remains un- 
altered. 

Correction of x axis orthogonality error is the final operation 
of the sample subroutine. Equation (3-1) is implemented at 
this point using a table lookup (for the sin function) , multi- 
plication and addition. The final result can be expressed as 

Hx = Hx' + Hy Sin e 


"^A byte is accepted terminology for an eight bit data quantity. 
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where the respective quantities have the following forms 

1 

Hy 
Sin e 

and 

Xi, X 2 are sign bits 

^■1 equal 0 or 1 depending on whether the respective 

i f J / is, 

term is to exist or not 

We can analyze the effects- of imprecision and truncation 
by noting that the sin function has eight significant binary 
bits resulting in a resolution of 1/256 or 90°/256 = 0,352°. 

The relative error in sin e is computed by Dahlquist 
[Ref. 4-12] as follows 

let a = the approximate value of sin e 
a = the exact value of sin e 

then the relative error in a is 

(a - a)/a if a / 0 

Since data in the sin table has been truncated, maximum 
relative error can be as large as +( 1 / 2 ^^) qj- +o.02%. 



From the definition of relative error we obtain the 
following relationships between exact, estimate and estimated 
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relative error 


a = a + ar = a(l + r) 

If ai, and 32 have relative errors of +0.39% and +0.02%, res- 
pectively, then 

3^32 = aid +0.0039) a (1 +0.00024) 

= aia 2 d +0.0039) (1 +0.00024) 

Thus, the relative error in a 133 is 

(1 +0.0039) (1 +0.00024) - 1 = 

+(0.0039) +(0.0039) (0.00024) +(0.00024) 

= +(0.0041) 

Since the maximum value of Sin e to be encountered occurs 
when the orthogonality error (s) is 1 degree, sin e - 0.017 
maximum. The maximum value for Hy can be 0.6 gauss or 2048 
units. Maximum error due to imprecision in the product HySins 
is then 

Er Max = (2048 x 0.017) (1 + .0041) - (2048 x 0.017) 

= 0.1427 units 

Since only the integer portion is retained in the final 
product, insignificant error can be attributed to imprecision 
of the sin e terra in this case. Orthogonality error will be 
adequate ly corrected . 

B) Subroutines ROTX and ROTY 

These subroutines were developed in Chapter III and 
implement the equation of 2-11 required to compute horizontal 
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X and y magnetic field components. Equations to be implemented 
by the respective subroutines are 

Hx = Hxm Cos (pitch) + Hym Sin (pitch) Sin (roll) 

+ Hzm Sin (pitch) Cos (roll) (4-12) 


and 


Hy = Hym Cos (roll) - Hzm Sin (roll) (4-13) 

where Hxm, Hym and Hzm are measured field components made avail- 
able from the magnetometer via the analog subsystem. 

Since the transcendental functions are implemented using, 
table lookup and are limited in precision to 8 bits, imprecision 
in these variables will dominate in generating error. In part- 
icular, the sin/c.os terms will have relative error in the order 
of +1/256 or +0.391 while the measured field data has relative 
uncertainty of only +1/4096 or +0.02%. Multiplications will 
result in addition of the bounds for the relative error as 
illustrated in section 4-4A above. 

The transcendental terms above are limited in magnitude 
to 1.0 maximum while the field measurements can be 0.60 gauss 
max. In this case the individual product terms of (4-12) and 
(4-13) can have maximum errors of 

Er = (2048) (1 + 0.0041) - (2048) = 8.4 units 

Errors in Hy and Hx (4-12 and 4-13) will be maximum when 
roll and pitch are at 45 degrees and the fields are equal. In 
this case the error in Hy will be 
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- . 0 ., 

EHy = 4.94 - 5.94 = 17 37 

J. 1 .87 units 

Hx - 17.8 units maximum 

It Should be noted tha^ +-h^ 

-a. at multipi.3 o. 45 


C) Subroutines COSY and SINY 


These two subroutines compute th« 
sensor (when projected onto the hof l^etween the x axis 

north-south horizontal vector of eartf ' 
first two opeations of these s.,K "ngnetic field. The 

Oision multiplication and divisilr^T^ 

involved are 12 bits in lencth a variables 

preserving 18 bits, no error is introJuceT"""''^ 

computes the dssiTeTTK^ZiTll^lJ'^ subroutines 

completing an associative table 

required is to match a given data ° The task 

- ~ r 
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ascertain the true address^. The function stored in tabular 
form is cos^O where 0 varies from 45 to 90 degrees. Maximum 
error will therefore be induced while attempting to locate 
solutions (angles near 90 degrees if inadequate precision is 
provided. Error in this. region due to resolution of tabular 
data can be examined by noting the entries in Table 4-8 


0 

Cos^0 

Most Significant 
Binary Bit (2~^) 

90 

0 


89 

0.000305 

12 

88 

0.001218 

9 

87 

0.00274 

8 


Table 4-8 Cos^0 AND MOST SIGNIFICANT BINARY DIGITS 

provided to indicate the relative magnitudes 'of Cos^O in the 
region of 0 = 90 degrees. We observe that the most signifi- 
cant binary digit affected at 89 degrees is binary decimal 
digit 12 implying that the resolution of or Hy^/Hh^ 

(the argument of Cos^0) must be accurate to at least 1/2^^ or 
0.024%. 

Considering the horizontal field of earth’s magnetic 
vector as observed in laboratory experimentation at this lat- 
itude, we note that Hh is 730 units. At a heading of 89 deg- 
rees, Hx = 7 30 Cos 89 == 12.7 units. The argument would there 
fore be 

ARG = Hx^/Hh^ = 0.000305 


^The procedure determines the relative address by linear 
interpolation, then selects the closest address as the 
required angle for the solution. 
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Since, the squaring and division operations are conducted 
in double precision, precision is preserved and the algorithm 
should be able to resolve heading to at least one degree over 
all portions of the compass - 

D) Errors Due to the Remaining Subroutines 

Since all of the remaining subroutines work with data that 
has been rounded to a precision representing 1 degree or better 
and the computations involve addition or subtraction in double 
precision binary or binary coded decimal (BCD) format, we note 
that there will be no further appreciable error due to trunca- 
tion or rounding. 

4-5 MEASUREMENT ERROR SUMMARY- 

Errors due to sensors and measurement of their respective 
outputs were discussed in sections 4-2 and 4-3 above. Before 
proceeding with the analysis of errors, the total signal in- 
accuracy due to contribution from the many sources above will 
be summarized in Table 4-10. Total instrument error can then 
be computed by considering the propagation and enhancement of 
these errors during the computation process. 

Since the errors in Table 4-9 are stochastically independ- 
ent, we can compute error for any given signal level by find- 
ing the RSS of respective error sources. In this manner, the 
instrument error can be evaluated by considering all input 
signals with errors superimposed to produce an erroneous com- 
putation of heading. 
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I 


PARAMETER 

ERROR 

COMMENT 

r 

1. Magnetometer 
Offset 

Orthogonality 

Noise 

Gain 

Linearity 

=^0 

Negligible 

+0.01% 

+0.01% 

Corrected by software 
Corrected by software 

Proport. to signal level 

1 

y 

2. Analog Subsystem 
Sampling 

Negligible 

Sampling rate & filter- 
ing adequate 

f 

\ 

Multiplexer 

Transfer Accuracy 
Settling Time 
Rate 

Input Leakage 

+0.01% 

==0 

Proport.to signal level 


Sample and Hold 
Acquisition 
Hold 

Aperature Delay 

Offset 

Gain 

Nonlinearity 

Feedthrough 

=^0 
4mV 
= 0 
= 0 
^0 

+0.01% 
+0 . 01% 

Corrected by software 
Corrected by software 
Proport.to signal level 
Proport.to signal level 
and frequency 

r 

f 

A/D Converter 
Resolution 
Accuracy 
Offset 
Gain 

+1.2mV 
+1.4mV 
= 0 
= 0 


\ 

\ 


Table 4-9 SENSOR AND ANALOG SUBSYSTEM ERROR SUMMARY 


r 
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4-6 SAMPLE ERROR ANALYSIS 


Orthogonality correction using the algorithmic method can 
be verified by computing expected error prior to correction and 
comparing measured system output with the error predicition. 
Assuming that the angle between the x and y sensors exceeds 90 
degrees as in Fig, 4-13, we can proceed to compute error by 
noting the following relationships 

Hx = Hh Cos 
riy = Hh Sin (ip) . 

Hxi = Hx Cos £ - Hy Sin s 

True Yaw = = Cos-^ 

Computed yaw 

HxCosg - HySine 

(HxCose - HySine) 2 + Hy^)^ 

HhCosijiCose - HhSiniJjSine 

[(HhCostj/Cose - HhSinti<Sine) 2-t. H^hSin^tl^]^ 

Using small angle approximations with e = 0.79° 

Cose = 1 and Sine 0,014 


= Cos 


= Cos 


-1 


-1 


then 


ijjm = Cos 


-1 


Hh Cost - 0,014 Sinii; 


[Hh^ (Cost - 0.014 Sint)^ + Hh^Sin^t^]^ 


Computed error 


Error = ti^ *" t 



100 


N 

X 


Y 

Fig. 4-13 ANGLE (X - Y)>90° 

We can now evaluate computed yaw angle given a part- 

icular yaw (i(j) and the horizontal field vector (Hh) . Heading 
error for horizontal field vector of 730 units at various yaw 
angles with pitch and roll angles of zero degrees is tabulated 
in Table 4-10 and plotted along with actual measured yaw error 
(data taken during experimentation of Chapter V) in Fig. 4-14. 


Heading 

(Degrees) 

Computed 

Error 

(Degrees) 

Heading 

(Degrees) 

Computed 
Error 
(Degrees ) 

90 

0.8 

290 

0.7 

70 

0.7 

270 

0.8 

50 

0.5 - 

250 

0.7 . 

30 

0.2 

230 

0.5 

10 

0.0 . 

210 

0.2 

350 

0.0 

190 

0.0 

330 

0.2 

170 

0.0 

310 

0.5 

150 

0,2 



130 

0.5 

Table 4-10 COMPUTED HEADING ERROR WITH Hh 

= 730 UNITS 




Coii'iputec 



Fig. 4-14 COMPARISON OF COMPUTED AND MEASURED ORTHOGONALITY INDUCED ERROR 
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4-7 CONCLUSIONS 

The preceeding error analysis has identified potential 
error sources along with relative magnitudes of error to be 
expected. Magnetometer sensor and analog subsystem errors 
were identified and analyzed individually. During this analy- 
sis it became apparent that errors due to sensor offset and 
nonorthogonality dominated and would severely limit total in- 
strument performance. The relative magnitudes of these errors 
and their mode of contribution would have degraded system cap- 
acity . 

By carefully characterizing the offset and orthogonality 
error it was determined that these systemmatic errors could be 
reduced by appropriate programming. A need to identify the 
extent of each error unique to the laboratory instrument im- 
posed a need to evaluate the instrument empirically. Using 
earth's magnetic field and the laboratory test fixture (de- 
scribed in Chapter V) to provide control inputs each of the 
parameters was identified and measured. An algorithm with the 
empirically determined correction coefficients was included in 
the final system to reduce the error and to improve final system 
performance. The remaining potential error sources were tab- 
ulated and relative magnitudes noted. 

Processing errors due to register precision and truncation 
were analyzed by considering pertinent subroutines individually. 
It was noted that the relative error bounds add when multiply- 
ing variables with relative error. In addition, it was noted 
that error accrued during processing is proportional to sensor 
signal levels involved. The final uncertainty is then propor- 
tional to actual aircraft attitude with error increasing as 
displacement from level flight occurs. Computational error is 
also noted to increase at particular headings causing the error 
function to peak at specific yaw angles. 
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A 

The sample error analysis clearly shows that a correlation 
between sensor nonorthogonality induced error and measured (un- 
corrected) data exists. By predicting and computing an error 
function prior to experimentally verifying the result we gain 
confidence that the sensor characteristics derived empirically 
in previous sections are correct. 
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CHAPTER V 

LABORATORY EVALUATION OF THE ATTITUDE INDEPENDENT 
REMOTE MAGNETIC INDICATOR AND HEADING INSTRUMENT 

5-1 INTRODUCTION 

• This chapter addresses laboratory evaluation of the micro- 
processor based computer designed to implement the heading 
measurement instrument. An integral part of this instrument 
was the three axis fluxgate magnetometer used to implement the 
attitude independent remote magnetic indicator of Chapter II. 

The laboratory evaluation was designed to investigate empiric- 
ally the effects of physical parameters that would otherwise 
be impossible to assess. 

Although phenomena such as noise/ magnetic field gradient/- 
sensor orthogonality errors and offset errors can be predicated, 
combined effects on the propo"sed instrument and remote magnetic 
indicator are best evaluated in the laboratory. In addition, 
it was noted that errors due to sensor offset and nonorthogon- 
ality could be corrected by software included with the sample 
subroutine. Determination of the effectiveness of this cor- 
rection technique necessitated laboratory measurements of the 
errors (to determine correction constants) and comparison of 
data prior to and following corrections. 

The chapter begins by discussing laboratory test apparatus 
designed to evaluate the instrument. Actual data measured and 
recorded during experimentation is then presented in both tab- 
ular and graphic form to facilitate comparison and evaluation. 
Finally, the laboratory data is discussed and it is concluded 
that the remote magnetic indicator used with the heading mea- 
surement instrument results in a viable alternative to conven- 
tional heading measurement systems. The microprocessor based 
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computer implentation of the instrument has added unique sensor 
measurement correction ability that enhances performance of 
otherwise marginal sensors. In this manner limitations in 
systems performance that now exist due to sensor inadequacy 
can be minimized without incurring the burden of using more 
expensive sensors . 

5-2 TEST APPARATUS 

A) Electronic Subsystem 

The microprocessor based computer (illustrated in photos 
5-1 and 5-2) was constructed on printed circuit boards consist- 
ing of a central processing card, two memory cards (2K bytes 
capacity each) and an output board. A separate analog subsystem 
card contained the multiplexer, ^sample and hold, analog to 
digital converter and trimming potentiometers. The circuit 
cards were all organized with edge connectors and mounted ver- 
tically into a hand wired backplane assembly as shown in photos 
5-1 and 5-2. 

The card in the left foreground of photo 5-1 served as the 
output display with three seven-segment displays displaying 
significant figures of system heading. A small printed circuit 
in the right foreground of photo 5-1 contained potentiometers 
used to generate analog signals proportional to roll and pitch 
signals (simulating gyroscope outputs). Cards shown vertically 
mounted in photo 5-2 can be identified from right to left as 
the analog subsystem, two memory cards and the central pro- 
cessing card. The large integrated circuit shown on the CPU 
card is the Signetics 2640 microprocessor. 

B) Sensor Assembly 


To evaluate the effects of combined aircraft pitch, roll 


REPRODUCIBILITY OF THE 
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Photo 5-1 MICROPROCESSOR BASED HEADING COMPUTER 



Photo 5-2 


CENTRAL PROCESSOR, MEMORY AND 
ANALOG SUBSYSTEM 
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and yaw a three axis gimbal apparatus was required. In add- 
ition, since angular measurements were required, a means of 
measuring angular rotation in each of the three exes was pro- 
vided. The gimbal apparatus as illustrated in photos 5-3 
and 5~4 was fitted with large protractors centered on 'the 
rotation axes. Pointers were provided to enable angular rot- 
ation measurements on the respective protractor scales. Since 
the angular precision on each protractor scale resolved angular 
position to 0.5 degrees, angular measurements to a resolution 
of at least 0.5 degrees were possible- Angular position was 
measured by estimating the decimal place of each measurement 
with accuracy to +0.5 degrees ensured. 

Since the three axis magnetometer (housed in the rectang- 
ular block of photos 5-3 and 5-4) measured ambient magnetic 
fields the test apparatus was constructed of nonferrous mat- 
erial. This ensured that local fields due to residual mag- 
netic fields in the test apparatus would be minimized. In 
addition, since the material had low permeability, there would 
be little deformation of the local field causing error due to 
changing field gradient. 

The sensor package shown in photos 5-3 and 5-4 was physic- 
ally mounted such that the sensors were centered as close to 
the center of the gimbal as possible. This precaution ensur- 
ed that measurement error due to sensor translation was minim- 
ized^ . During instrument evaluation, the entire gimbal assembly 
and sensor were leveled and mounted in a Helmholtz coil assembly 
as illustrated in photo 5-5. Although the coils were not acti- 
vated during experimentation, the rotations in heading were 

^Since the local magnetic field has a nonzero gradient, field 
measurements include a component due to translation of the 
sensor axes. This component of measurement produces unaccept- 
able error in a system designed to measure field components 
that change due to rotation. 
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Photo 5-3 


MAGNETOMETER SENSOR MOUNTED ON 
GIMBALLED TEST FIXTURE 



Photo 5-4 SENSOR AND 


gimbae assembly with protractors 



original ?A 
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carefully controlled since the gimbal assembly was an integral 
part of the Helmholtz coil fixture with the vertical rotation 
axis serving as the system yaw axis. 



Photo 5-5 TEST FIXTURE MOUNTED IN HELMHOLTZ COIL ASSEMBLY 


5-3 HEADING MEASUREMENTS WITH NO OFFSET CORRECTION 

By- maintaining heading of the text fixture constant (no 
rotation about the vertical axis) and varying both pitch and 
roll angle, the instrument display was observed to vary. This 
variation gave a direct measure of instrument error since a 
constant heading was maintained arid a constant display was to 
be expected. 

Data variations were recorded in Tables 5-1 and 5-2 and 
plotted on Figures 5-1 and 5-2. With only +10 degree variation 
in pitch combined with +30 degree variation in roll we note 
that the heading display varies 14 degrees. Obviously, 


Roll Angle 
(Degrees) 


0 



Table 5-1 


HEADING COMPUTED AT A FIXED YAW ANGLE 
WITH NO OFFSET CORRECTION 



Table 5-2 


heading computed at a fixed yaw angle 

WITH NO OFFSET CORRECTION 
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instrument operation indicated excessive error requiring more 
elaborate sensors or correction of a sensor inadequacy. 

5-4 HEADING MEASUREMENTS TO INVESTIGATE ORTHOGONALITY ERROR 


System performance was evaluated by initially aligning the 
sensors with zero pitch and roll angle. Sensor Z was position- 
ed vertically with positive direction downwards. By observing 
the Z axis output^ as the test fixture was rotated about the 
vertical axis, adjustments were made in pitch and roll angle 
to minimize coning of the Z axis. Angular measurements on 
the respective roll and pitch axis protractors were then made 
to establish the initial reference attitude angles. 

Heading measurement accuracy was evaluated by rotating 
the test fixture in the horizontal plane until the display 
flickered between (XX9) and (XX9+1) . The rotation was then 
continued .a very small amount until a steady display (multiple 
of 10 degrees) was observed^. Measurements ranging- from 0 to 
350 degrees were made by recording angular position required 
to produce specific heading data displays. Sets of data were 
recorded at various combinations of pitch and roll then tab- 
ulated in Tables 5-3 through 5-8. Relative error was computed 
by determining angular position expected at each display value 
and then computing the difference in angular positions. Errors 
at the roll extremes -of +44 degrees are plotted for pitch 
angles of plus and minus 20 degrees on Fig- 5-3 through 5-6 
inclusive. 


^A special subroutine was. used to display Z axis data directly 
in BCD format on the seven bar output display. 

^This measurement technique ensured that all heading measure- 
ments were made identically. In addition, error due to, system 
imprecision was reduced. 
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Data in Tables 5-5, 5-6 and Fig. 5-3, 5-4 were recorded 
with no sensor orthogonality error correction implemented. Data 
in Tables 5-7, 5-8 and Fig. 5-5, 5-6 was recorded with the 
sensor orthogonality correction implemented. Comparison of 
these data indicate that considerable improvement in accuracy 
is achieved by correcting sensor orthogonality error. 

5-5 CONCLUSIONS 

Laboratory evaluation of the heading measurement instru- 
ment has shown that the algorithms developed in previous 
chapters are viable. Operation of the device in a laboratory 
environment has enabled empirical evaluation of the system 
under adverse combinations of noise, field gradient and sensor 
plus instrument error sources. 

Test apparatus described in section 5-2 served to. enable 
controlled simulation of roll, pitch and yaw rotations. The 
apparatus was nonmagnetic in nature and contributed insignifi- 
cant error due to field pertebation. Mounting of protractors 
and pointers on the test apparatus made angular measurements 
possible to a precision of at least +0.5 degrees. 

Effects of sensor offsets were evaluated in section 5-3 by 
recording system heading computations when only roll and pitch 
varied. Since the variations in Figures 5-1 and 5-2 prior to 
offset correction exceed the maximum excursions of Figures 5-3 
and 5-4 by at least a factor of two (angular excursions in 
first set also less than in the record) and we note that offset 
errors were corrected prior to recording data in the second set 
of data, we conclude that offset in magnetometers can be a 
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PITCH ANGLE 0 DEGREES 
ROLL ANGLE 0 DEGREES 


Heading 

Displayed 

(Degrees) 

Angular 

Position 

(Degrees) 

Relative 

Error 

(Degrees) 

10 

275.3 

0.3 

30 

295.2 

0.2 

50 

315.4 

0.4 

70 

335.5 

0.5 

90 

355.5 

0.5 

130 

34.8 

-0.2 

150 

54.3 

-0.7 

170 

74.0 

-1.0 

190 

94.0 

-1.0 

210 

114.4 

-0.6 

230 

■135.0 

0.0 

250 

155.5 

0.5 

270 

174.7 

-0.3 

290 

195.5 

0.5 

310 

215.2 

0.2 

330 

235.2 

0.2 

350 

254.7 

-0.3 


Table 5-3 


REFERENCE DATA MEASUREMENTS OF HEADING 
TAKEN WITH NO ORTHOGONALITY CORRECTION 



PITCH ANGLE 0 DEGREES 
ROLL ANGLE 0 DEGREES 


Heading 

Angular 

Relative 

Displayed 

Position 

Error 

(Degrees ) 

(Degrees) 

(Degrees) 

10 

276.6 

1 

o 

30 

296.6 

-0.4 

40 

306.7 

-0.3 

50 

316.7 

-0.3 

60 

327.2 

+ 0.2 

70 

337.0 

0.0 

90 

355.9 

-1.1 

130 

37.0 

0.0 

160 

67.2 

0.2 

190 

96,8 

-0.2 

220 

126.9 

-0.1 

250 

157.0 

0.0 

290 

186,9 

-0.1 

310 

216,6 

-0.4 

340 

246.9 

-0.1 

350 

256.4 

-0.6 


Table 5-4 HEADING MEASUREMENTS WITH OFFSET AND 
ORTHOGONALITY CORRECTIONS MADE 



Roll = 44° Roll = 20° Roll = -20° Roll = -44 
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TABLE 5-5 HEADING MEASUREMENTS AT PITCH = 20° WITH NO 
ORTHOGONALITY CORRECTION 
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Heading 

Displayed 

(Degrees) 

20 
40 
60 
80 
90 
140 
160 
180 
200 
220 
240 
260 
280 
300 
320 
340 
0 


Roll = 

0° 

Roll = 

Angular 

Position 

Error 

Angular 

Position 

287.3 

0.3 

287.5 

307.2 

0.2 

307.5 

327.0 

0.0 

327.9 

346.3 

-0.7 

347.5 

355.9 

-1.1 

357.0 

46.9 

-0.1 

47.0 

67.0 

0.0 

67.3 

86.6 

-0.4 

86.5 

107.0 

0.0 

107.3 

127.6 

0 . 6 

127.3 

147.2 

0.2 

147.4 

167.0 

0.0 

167.1 

187.3 

0.3 

187.1 

207.2 

0.2 

207.1 

227.3 

0.3 

227.0 

247.4 

0.4 

247.3 

266.4 

-0.6 

267.0 


Roll = '44 


irror 

Angular 

Position 

Error 

0.3 

286.3 

-0.7 

0.3 

306.0 

-1.0 

0.7 

325.9 . 

-1.1 

0.3 

345.6 

-1.4 

0.2 

356.2 

-1.8 

-0.2 

46 . 5 

-0.5 

0.1 

67.1 

0.1 

-0.7 

87.0 

0-0 

0.1 

108.1 

1.1 

0.1 

128.3 

t — 1 

0.2 

148.5 

1.5 

-0.1 

168.5 

1.5 

-0.1 

188.5 

1.5 

-0.1 

208.0 

O 
1 — 1 

-0.2 

227.5 

0.5 

• 

o 

+ 

247.4 

0.4 

-0.2 

266.0 

-l.C 


Table 5-7 


miNG measurements AT 

AMH ORTHOGONALITY 


MADE 


= 20 ° 

CORRECTION 
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Heading 

Displayed 

(Degrees) 

Roll = 44° 
Angular 

Position Error 

Roll = - 

Angular 

Position 

-44° 

Error 

20 

287.0 

-0.4 

286.5 

-0.3 

40 

307.0 

-0.4 

306.7 

-0.1 

60 

327.0 

-0.4 

326.8 

0.0 

80 

347.0 

-0.4 

346.8 

0.0 

90 

356.7 

-0.7 

355.6 

-0.2 

140 

47.8 

0.4 ■ 

48.0 

1.2 

160 

67.7 

0.3 

68.0 

1.2 

180 

87.0 

-0.4 

87.8 

1.0 

200 

107.6 

0.2 

108.0 

1.2 

220 

127.7 

0.3 

128.0 

1.2 

240 

147.8 

0.4 

147.4 

0.6 

.. 260 

167.6 

0.2 

167.2 

0.4 

280 

137.6 

0.2 

187.0 ■ 

0.2 

300 

207.1 

-0.3 

206.5 

-0.3 

320 

227.1 

-0.3 

226.6 

-0.2 

340 

247.1 

-0.3 

246.9 

0.1 

0 

266.3 

-0.9 

266.0 

-0.8 


TABLE 5-8 


HEADING MEASUREMENTS AT PITCH = -20° WITH 
OFFSET AND ORTHOGONALITY ERROR CORRECTED 
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major error source‘s. Additionally, we note that the .correction 
of off-set error in sensors has been successful. Experimental 
results have verified that not only can offset errors be deter- 
minded (Chapter IV) , but a suitable algorithm can be implement- 
ed in the computer to improve system operation. It is postu- 
lated that offset error correction can be extended to include 
correction of varying offset values (functions of temperature 
and supply voltage) by monitoring error causing variables 
(example temperature) and computing correction constants prior 
to offset correction as above. 

Errors induced by sensor nonorthogonality were predicted 
in Chapter IVsection 4-2 and verified by plotting expected 
error along with measured error in Fig. 4-14. The curves of 
Fig. 4-14 were plotted for heading rotations with no pitch or 
roll angle. To evaluate' system performance and the effect of 
orthogonality error with combined angular rotations, measure- 
ments of heading error were plotted in Fig. 5-3 thorugh 5-6 
inclusive. 

Comparison of these data indicate that maximum excursions 
of error as a function of heading are significantly less when 
orthogonality corrections are made. It is also postulated that 
data could be improved further by similarly correcting ortho- 
gonality error in te 2 axis sensor^. 

In summary, the experimental evaluation has provided in- 
sight into the operation of an attitude independent remote mag- 
netic indicator and heading computer in the '‘real world" 


‘^This corroborates the observations predicted during error 
analysis in Chapter IV. 

^We note that the error excursions are 'functions of pitch 
and roll and that Z axis data is used in the rotation algorithm. 
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environment complete with all contributing error sources. The 
error analysis evolved during development of the system has 
proven adequate in that an operational system was developed. 
Major error sources were measurable as predicted and the means 
of reducing their effects were successfully implemented. Cor- 
rection of sensor offset and orthogonality error required an 
empirical evaluation of the respective sensor. These evalua- 
tions were performed, the errors characterized, correction co- 
efficients determined, and correction algor itms implemented. 

Successful implementation of these corrections was evi- 
denced by significant reductions in system error. The correct- 
ion methods presented can be extended in future with the net 
result that less demand is required of physical sensors if the 
sensor parameters can be established empirically prior to 
completion of instrument design. Utilization of a microproces- 
sor in the instrument has added the computational flexibility 
required to facilitate accommodation of sensors with varying 
error magnitudes. 


appendix a 


lists the instruction set of the Signetics 
Thrs appendix 1 implement the heading instru- 

2650 microprocessor chip 

ment. 
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APPENDIX B 


This appendix contains alisting of the assembly language 
program used to implement the remote magnetic indicator head- 
ing algorithm. The program was assembled on ^ the A2650 cross 
assembler program operational on the HP 2100 ' computer at the 
University of Santa Clara. 
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APPENDIX C 


The transcendental functions used throughout the heading 
computation algorithm were implemented using a table look up 
procedure. To generate the respective look up tables in com- 
puter memory data was first generated using algol programs. 
This technique expedited modifications to tabular data and 
provided output data in a convenient (hexadecimal) format. 

Programs that calculated Cos (0.) and Cos^ (0)to eight 
bit and sixteen bit resolution respectively are included. 
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